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ABSTRACT 


The  theoretical  study  presented  is  concerned  with  predicting  the  flow 
details  in  the  region  immediately  behind  a  slender  wedge  moving  at  supersonic 
speed.  An  asymptotic  solution  is  obtained  for  the  laminar  near  wake,  which  is 
valid  in  the  limit  as  Reynolds  number  based  on  external  flow  properties  tends 
to  infinity,  with  the  external  Mach  number  held  fixed.  In  this  limit,  the  interi¬ 
or  of  the  recirculating  flow  is  approximately  incompressible  and  inviscid;  the 
viscous  layer  is  a  continuous  thin  layer  with  large  transverse  velocity  gradient; 
and  the  typical  velocity  of  the  recirculating  flow  must  be  of  the  same  order  as 
that  of  the  external  flow,  although  numerically  it  is  small . 

For  the  inviscid  core,  it  is  shown  from  consideration  of  the  dissipative 
terms  that  the  temperature  and  vorticity  are  nearly  constant.  These  results 
demonstrate  the  fact  that,  even  though  the  core  flow  i  inviscid  and  nonconduct¬ 
ing,  and  the  primary  effects  of  diffusive  exchanges  are  limited  to  the  continuous 
viscous  layer  enclosing  the  recirculating  flow,  the  accumulative  effects  of  very 
small  diffusive  exchanges  are  experienced  throughout  the  entire  recirculation 
region. 


There  exists  an  equilibrium  vorticity  for  the  inviscid  core  so  as  to  satisfy 
the  condition  of  conservation  of  angular  momentum  for  the  recirculating  flow.  At 
the  rear  stagnation  point,  a  tentative  reattachment  condition  similar  to  Chapman’s 
is  used  but  based  on  the  minimum  pressure  increase  attainable  at  the  rear  stag¬ 
nation  ooint  from  the  turning  of  a  streamline  in  the  shear  layer  to  a  direction  par¬ 
allel  to  the  wake  center  line.  Based  on  these  two  conditions,  a  unique  solution 
for  the  near  wake  flow  can  be  obtained  for  a  given  upstream  flow  condition.  Nu¬ 
merical  results  are  in  general  agreement  with  the  available  experimental  data . 
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THE  LAMINAR  NEAR  WAKE  BEHIND  A  SLENDER  WEDGE 
AT  SUPERSONIC  SPEED  AND  HIGH  REYNOLDS  NUMBER* 

I 

INTRODUCTION 

Renewed  effort  has  been  directed  recently  to  the  theoretical  and  experi 
mental  investigation  of  the  wake  behind  a  body  moving  at  high  speeds  through 
air.  Considerable  progress  has  been  made  in  analytical  studies  of  both  the 
laminar  and  turbulent  far  wake  (e.g.,  see  the  review  paper  of  Lykoudis1)-  How 
ever,  accurate  prediction  of  the  far-wake  flow  properties  requires  knowledge 
of  these  same  properties  at  a  distance  of  a  few  body  heights  behind  the  body. 
Because  of  the  extreme  complexity  of  the  flow  in  the  near -wake  region,  the 
attempts  to  obtain  satisfactory  analytical  results  have  been  only  partially  suc¬ 
cessful.  Accordingly  current  studies  are  focusing  on  the  near  wake  in  the  hope 
that  an  adequate  analysis  of  the  complex  coupling  processes  occurring  within 
this  region  can  be  carried  out.  The  present  investigation  is  concerned  with  the 
laminar  near  wake  of  a  slender  wedge  at  low  supersonic  Mach  numbers  and  high 
Reynolds  numbers. 

Knowledge  of  experimental  results  is  essential  to  the  developing  of  a 
mathematical  model,  both  as  a  guide  for  the  choice  of  approximations  and  for 
a  check  on  the  accuracy  of  the  results  predicted  using  the  model.  At  the  base 
of  the  wedge  it  is  known  that  the  external  flow  expands  to  a  lower  pressure, 
through  a  Prandtl -Meyer  expansion.  The  boundary  layers  along  the  wedge  sur¬ 
faces  separate  and  form  free  shear  layers  which  enclose  a  region  of  steady  re¬ 
circulating  flow.  These  results  are  evident,  for  example,  in  the  experiments 

*  This  report  is  based  on  a  dissertation  submitted  in  partial  fulfillment 

of  the  requirements  for  the  degree  of  Doctor  of  Philosophy  at  The  University 
of  Michigan,  1968. 
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carried  out  by  Hama  ,  for  various  Mach  numbers  between  about  2.0  and  4.0, 
and  for  a  range  of  Reynolds  numbers  such  that  the  free  shear  layer  was  lami¬ 
nar  at  t he  lower  values.  Among  the  results,  it  was  shown  that  the  pressure 
is  almost  constant  on  the  base;  the  dividing  streamline  is  nearly  straight  over 
most  of  its  length  from  the  base  to  the  rear  stagnation  point}  and  a  "lip  shock" 
originates  in  the  shear  layer,  very  close  to  the  comer.  This  lip  shock  may 
have  considerable  strength  for  flows  with  a  large  expansion  ratio  at  the  corner, 

but  remains  weak  for  flow  over  a  slender  wedge  at  low  super  sonic  Mach  num- 

3  4  o 

bers  .  Batt  has  made  measurements  in  the  wake  behind  a  20  wedge  atMo  =  6 

and  at  various  Reynolds  numbers  such  that  the  free  shear  layer  is  laminar. 
The  results  indicate  that  the  location  of  the  rear  stagnation  point  is  not  very 
sensitive  to  the  variation  in  Reynolds  number,  and  that  the  recirculating  flow 
has  a  very  low  density.  Martellucci  et  4l  have  studied  the  turbulent  near 
wake  of  a  cone  at  Ma  =  6.  Their  results  show  that  the  pressure  at  the  base  is 
nearly  constant, but  the  velocity  at  the  axis  in  the  recirculation  region  may  cor¬ 
respond  to  Mach  numbers  as  high  as  0.8 .  For  a  Mach  number  equal  to  16 
6 

Todisco  and  Pallone  found  that  a  region  of  constant  temperature  exists  in  the 
recirculating  flow  of  the  laminar  near  wake  behind  a  slender  wedge  or  cone. 

These  and  other  experimental  results  therefore  suggest  that  the  near 
wake  of  a  slender  body  can  be  divided  into  regions  of  different  physical  charac¬ 
ter.  An  external  inviscid  flow  results  from  the  comer  expansion  of  the  flow 
outside  the  boundary  layer.  A  free  shear  layer  is  formed  by  the  expanded 
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boundary  layer  which  mixes  wit*  the  relatively  low-speed  recirculating  flow. 
The  velocity  gradient  is  large  nttr  the  dividing  streamline,  which  separates 
the  recirculating  fluid  from  the  fluid  which  passed  over  the  wedge.  Since  the 
flow  near  the  dividing  streamline  is  directed  toward  the  wake  center  line,  an 
increase  in  pressure  occurs  in  the  vicinity  of  the  rear  stagnation  point,  in  or¬ 
der  that  the  flow  may  be  turned.  The  fluid  with  total  pressure  higher  than  on 
the  dividing  streamline  escapes  downstream  and  forms  the  far  wake,  whereas 
the  fluid  with  lower  total  pressure  turns  back  toward  the  base  and  constitutes 
the  recirculating  flow.  The  shear  stress  along  the  dividing  streamline  drives 
the  recirculating  flow  against  the  retarding  force  along  the  base  of  the  wedge . 

Since  the  near -wake  flow  field  is  determined  by  a  complicated  inter¬ 
action  between  the  flow  in  the  shear  layers,  the  inviscid  external  flow,  and  the 
internal  recirculating  flow,  analytical  descriptions  have  only  been  obtained  by 
some  rather  rough  approximation  or  through  a  semi -empirical  formulation. 

Several  simplified  near -wake  models  have  been  proposed. 

7  8  9 

The  Chapman  ’  ’  dividing-streamline  model  is  based  on  the  conserve 
tion  of  mass  within  the  recirculating-flow  region .  The  mass  flow  in  the  shear 
layer  below  the  dividing  streamline  is  balanced  by  the  mass  flow  reversed  by 
the  pressure  rise  through  the  recompression  zone  near  the  rear  stagnation 
point.  The  base  pressure  is  determined  by  requiring  that  the  total  pressure 
along  the  dividing  streamline  be  equal  to  the  final  static  pressure,  and  the  total 
pressure  is  assumed  to  be  fully  recovered  at  the  rear  stagnation  poin'.  It  is 
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assumed  that  the  viscous  shear  layer  has  zero  initial  thickness.  The  external 
flow  mixes  with  the  recirculating  flow,  which  is  considered  to  have  zero  veloc¬ 
ity,  and  it  is  assumed  that  the  thickness  of  the  viscous  shear  layer  remains 
small  in  comparison  with  the  dimensions  of  the  recirculation  region.  It  is  fur¬ 
ther  assumed  that  the  dividing  streamline  is  straight  and  that  the  boundary- 
layer  equations  may  be  employed  to  solve  the  shear -layer  problem.  Baum, 
King  and  Denison1*^  improve  on  Chapman's  model  by  considering  a  nonzero  ini¬ 
tial  viscous -layer  thickness-iT  They  propose  that  the  separated  shear  layer  is 
a  thin  mixing  region  between  the  outer  flow  and  the  essentially  stagnant  recir¬ 
culating  flow  with  the  suction  of  fluid  along  the  inner  edge  of  the  shear  layer 
made  up  by  a  source  introduced  at  the  rear  stagnation  point.  Taba^  makes  a 
further  improvement  by  taking  into  consideration  not  only  the  finite  initial  vis  - 
cous -layer  thickness  but  also  the  trailing-edge  expansion.  The  shear  layer  is 
regarded  as  consisting  of  an  outer  shear  layer  where  inertia  effects  are  domi¬ 
nant  and  an  inner  mixing  layer  where  both  viscous  and  inertia  forces  are  im¬ 
portant.  The  velocity  at  the  inner  edge  of  the  shear  layer  is  again  assumed 
zero . 

12 

A  different  model  is  used  by  Reeves  and  Lees  to  analyze  the  near 

wake  of  a  circular  cylinder .  The  boundary -layer  assumptions  are  regarded  as 

13 

valid  for  the  flow  in  the  near  wake.  Stewartson's  reverse -flow  profiles,  ob¬ 
tained  by  an  extension  of  the  Falkner-Skan  similar  solutions,  are  adopted  to 
describe  the  region  with  reversed  flow.  In  this  application  the  profiles  are 
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characterized  by  two  length  scales,  the  base  height  and  the  boundary -layer 

thickness  at  separation.  A  two-moment  integral  method  is  used  to  obtain  the 

14 

flow  properties .  Golik,  Webb,  and  Lees  extend  the  method  to  study  the  near 
wake  behind  a  wedge.  A  three -moment  integral  method  is  employed,  in  order 
that  i  ne  more  degree  of  freedom  may  be  available  in  representing  the  flow,  and 
the  axis  velocity  is  uncoupled  from  the  Stewartson  profiles .  It  is  discovered 
that  only  one  value  for  die  base  pressure  allows  the  integral  curve  of  the  solu¬ 
tion  to  pass  through  a  saddle  point  of  the  system .  Hie  requirement  that  the 
solution  curve  pass  through  this  critical  point  is  regarded  as  a  uniqueness  cri  - 
terion  for  the  near  wake. 

Dean1"*  and  Oswatitsch^  employ  a  double  power-  feries  expansion  of  the 

stream  function  about  the  separation  and  reattachment  points  for  the  study  of  the 

local  behavior  of  incompressible  and  viscous  flows  near  these  points .  These 

1 7  18 

local  solutions  are  extended  by  Reeves  and  Buss  ’  to  the  study  of  the  near 
wake  of  a  slender  wedge,  which  is  regarded  as  controlled  by  a  thin  shear  layer 
of  relatively  low  energy  fluid.  The  stream  function  is  expanded  about  the  rear 
stagnation  point  and  the  coefficients  of  the  series  can  be  determined  by  substi¬ 
tuting  in  the  full  Navier -Stokes  equations,  satisfying  suitable  conditions  at  sep¬ 
aration  and  at  the  rear  stagnation  point,  and  matching  with  die  downstream  flow . 
Accordingly  it  is  found  that  the  geometrical  and  dynamical  properties  of  the  near 
wake  are  determined  by  the  local  Reynolds  number  and  the  velocity  profile  at  the 


rear  stagnation  point. 
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Weiss  suggests  a  model  which  consists  of  three  coupled  regions:  a 
vortical  outer  flow,  which  results  from  an  approximately  inviscid  expansion 
of  the  boundary  layer  at  the  trailing  edge  of  the  wedge,  is  treated  by  the  meth¬ 
od  of  characteristics;  the  flow  in  a  thin  viscous  layer  above  the  dividing 
streamline,  which  has  large  velocity  gradient  and  can  be  described  by  the 
boundary -layer  equations,  is  solved  by  a  modified  Oseen  linearization  tech¬ 
nique;  and  the  recirculating  flow,  which  contains  flow  of  relatively  low  veloc¬ 
ity  and  is  described  by  the  full  Navier -Stokes  equations,  is  solved  by  a  finite - 
difference  technique.  A  numerical  iterative  scheme  is  proposed  such  that  the 
flow  properties  at  the  boundaries  of  the  coupled  region  are  properly  matched . 
These  conditions  are  obtained  by  demanding  that  velocity,  shear,  pressure, 

temperature  and  heat  transfer  be  continuous  for  any  given  base  pressure. 

20  21 

Prandtl  ’  suggests  that,  in  steady  incompressible  laminar  flows  with 
closed  streamlines  at  large  Reynolds  number,  viscous  effects  can  be  disre¬ 
garded  wherever  there  are  no  large  velocity  gradients  or  accumulative  effects  . 
Thus  the  viscous  effects  are  limited  to  narrow  regions  in  the  vicinity  of  bound¬ 
ary  surfaces  which  enclose  the  recirculating  flow.  Since  the  net  work  done  on 

a  fluid  particle  rgoving  along  any  closed  streamline  must  be  zero  for  a  complete 
. .  .  * 

cycle,  the  vorticity  can  be  shown  to  be  distributed  uniformly  in  the  interior  of 

the  recirculation  region  of  a  two  dimensional  flow  where  convective  effects  are 
22  23 

dominant.  Batchelor  ’  obtains  this  result  in  essentially  the  same  way,  by 
taking  the  line  integral  of  the  complete  equation  of  motion  around  a  closed 
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streamline  He  pro,.>oses  that  the  value  of  vorticity  may  be  determined  by  the 
condition  that  the  boundary  layers  enclosing  this  recirculating  flow  must  also 
be  in  steady  motion.  This  implies  that  in  incompressible  and  inviscid  wakes 
behind  two-dimensional  bodies  the  vorticity  must  take  such  a  value  that  the 
relative  momentum  of  the  viscous  layer  enclosing  the  recirculating  flow  is  in¬ 
creased  over  part  of  its  path  and  decreased  over  other  parts,  the  net  effect  be¬ 
ing  to  allow  the  viscous  layer  to  have  a  steady  motion.  It  is  implied  that  as 
Reynolds  number  becomes  large  the  velocity  in  the  recirculation  region  re¬ 
mains  of  the  same  order  as  the  velocity  in  the  external  flow. 

24 

Lenard  applies  Batchelor  s  idea  to  the  near  wake  of  a  slender  cone  at 
large  Reynolds  number  and  supersonic  Mach  number.  The  near  wake  is  divid¬ 
ed  into  the  inviscid  core  region,  the  corner -expansion  region,  the  external  in- 
viscid-flow  region,  the  recompression  region  and  tlm  viscous-layer  region, 
which  in  turn  is  subdivided  into  the  base  boundary  layer  and  the  free  shear 
layer.  The  base  boundary  layer  on  the  base  surface  of  the  cone  is  represented 
by  an  axially  symmetric  stagnation  point  boundary  layer  with  highly  vortical 
outer  flow.  The  free  shear  layer  which  grows  into  the  inviscid  external  flow 
and  inner  recirculating  flow  is  described  by  the  von  Mises  form  of  the  bound¬ 
ary  layer  equations .  The  external  flow  i«  considered  as  an  isentropic  expan¬ 
sion  of  the  flow  upstream  of  the  corner  to  the  base  pressure.  The  inviscid  in¬ 
ner  recirculating  flow  is  regarded  as  having  constant  properties  along  stream¬ 
lines.  The  integral  relations  which  are  obtained  by  integration  of  the  govern- 
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ing  differential  equations  along  a  closed  streamline  lying  entirely  within  the 
viscous  layers  provide  the  necessary  conditions  for  a  unique  solution.  A  nu¬ 
merical  procedure  is  suggested  for  obtaining  a  solution  to  the  equations . 

25 

Guderiey  and  Greene  ‘  adopt  Batchelor’s  proposal  for  the  study  of 
steady  compressible  swirl  flow  with  closed  streamlines  at  high  Reynolds  num¬ 
ber.  Some  special  conditions,  which  are  determined  by  the  dissipative  terms 
and  are  deduced  from  the  mathematical  conditions  for  solvability  of  the  equa¬ 
tions  for  approximations  of  a  higher  order,  can  be  obtained  for  the  inviscid 
core  flow.  These  conditions  have  the  physical  meaning  of  conservation  of  mo¬ 
mentum,  energy,  and  entropy  for  any  control  surface  generated  by  rotating  a 

streamline  in  the  recirculating  flow  around  the  axis  of  symmetry. 

14 

The  methods  proposed  by  Golik,  Webb,  and  Lees,  by  Reeves  and 
17  18  19 

Buss,  ’  and  by  Weiss  have  all  provided  results  which  compared  favorably 
with  experiment  for  one  or  more  combinations  of  Mach  number  and  Reynolds 
number.  However,  a  difficulty  remains  ineachcase  because  there  is  no  obvi¬ 
ous  procedure  for  estimating  the  error  or  for  systematically  improving  the 
approximation .  Therefore  it  appears  of  interest  also  to  study  asymptotic 
approximations  for  large  Reynolds  number  which  in  principle  can  be  extended 
to  include  higher -order  terms.  If  a  complete  solution  could  be  obtained  in  this 
manner,  it  is  expected  that  the  flow  details  would  be  predicted  more  accurately 
than  by  the  solutions  obtained  in  a  more  intuitive  way.  The  range  of  validity  of 
an  asymptotic  solution  would  probably  be  smaller,  but  the  solution  could  pro- 
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vide  a  useful  standard  of  comparison  for  solutions  obtained  by  the  other  meth¬ 
ods  quoted. 

23 

In  the  present  study  the  idea  proposed  by  Batchelor  for  imcompress- 

24  25 

ible  wakes  is  adopted,  in  the  manner  of  Lenard  and  Guderiey  and  Greene, 
to  provide  a  model  for  the  near  wake  behind  a  slender  wedge  at  supersonic 
speed.  That  is,  an  attempt  is  made  to  obtain  an  asymptotic  solution  for  the 
laminar  near  wake,  which  is  valid  in  the  limit  as  the  Reynolds  number  based 
on  external  flow  properties  tends  to  infinity,  with  the  external  Mach  number 
held  fixed.  Terms  considered  to  be  of  higher  order  are  not  obtained  here,  but 
it  is  pointed  out  that  a  complete  understanding  of  the  problem  requires  know¬ 
ledge  of  additional  terms  in  the  asymptotic  expansions.  It  is  important  to  em¬ 
phasize  that  a  solution  obtained  in  diis  limit  is  not  expected  to  describe  effects 
associated  with  high  Mach  number.  The  solutions,  of  course,  are  to  be  used 
not  at  arbitrarily  large  Reynolds  numbers  but  at  Reynolds  numbers  such  that 
the  flow  actually  is  laminar. 

As  suggested  by  Batchelor,  it  is  expected  that  as  the  Reynolds  number 
becomes  large  the  velocity  in  the  recirculation  region  remains  of  the  same 
order  as  the  velocity  in  the  external  flow,  although  numerically  smaller.  In 
this  limit,  therefore,  the  Reynolds  number  based  on  typical  properties  for  the 
recirculating  flow  will  also  approach  infinity,  although  in  any  application  the 
numerical  value  will  be  far  smaller  than  the  Reynolds  number  based  on  exter¬ 
nal  flow  properties.  It  is  believed  nevertheless  that  the  results  will  provide  a 
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useful  approximation  at  the  larger  values  of  Reynolds  number  for  which  the  flow 
remains  laminar . 

Since  experiments  suggest  that  the  Mach  number  in  the  recirculation 
region  may  not  be  very  large,  the  recirculation  region  is  assumed  to  have  con¬ 
stant  density,  and  an  estimate  is  made  to  show  that  the  effects  of  nonzero  Mach 
number  may  reasonably  be  neglected.  It  is  therefore  also  assumed  that  the 
dividing  streamline  is  straight  and  that  the  pressure  is  nearly  constant  except 
in  regions  near  the  corner  and  near  the  rear  stagnation  point.  In  the  limit  it 
is  assumed  that  the  dimensions  of  these  regions  approach  zero.  Since  the  size 
of  these  regions  will  increase  with  increasing  Mach  number,  it  seems  clear 
that  the  hypersonic  problem  would  have  to  be  studied  by  some  different  limit 
process  such  that  the  Mac  number  also  approaches  infinity  and  the  size  of 
these  regions  does  not  approach  zero.  Finally,  for  simplicity  zero  heat  trans¬ 
fer  to  the  wedge  surfaces  is  assumed,  and  the  Prandtl  number  is  taken  equal  to 
one. 

Solutions  are  derived  for  the  recirculating  flow,  for  a  shear  sublayer 
and  for  the  outer  part  of  the  free  shear  layer,  and  for  the  base  boundary  layer. 
Two  constants  remain  undetermined  in  these  solutions,  namely  the  base  pres¬ 
sure  and  the  vorticity  of  the  recirculating  flow.  Therefore  two  additional  con¬ 
ditions  are  needed.  A  balance  is  imposed  between  the  torques  exerted  on  the 
recirculating  flow  by  the  shear  stress  along  the  dividing  streamline  and  the 
shear  stress  along  the  base.  At  the  rear  stagnation  point  a  criterion  similar 
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CO  Chapman's  is  used,  but  based  on  flow  conditions  for  a  particular  streamline  inside 
the  shear  layer,  rather  than  on  conditions  outside  me  shear  layer.  It  is  pointed 
out  that  these  conditions  are  considered  as  tentative,  since  they  are  not  derived 
from  the  solutions  obtained.  If  second  approximations  could  be  obtained  in  the 
asymptotic  expansions  of  the  flow  properties  for  large  Reynolds  number,  it  is  be¬ 
lieved  that  these  or  similar  criteria  could  be  proposed  with  greater  confidence. 
Numerical  results  are  given  for  some  of  the  cases  studied  by  Hama,  and  rather 
good  agreement  is  found  between  predicted  and  measured  values  of  the  base 


pressure. 


II 


DESCRIPTION  OF  THE  PHYSICAL  PROBLEM 
AND  PROPOSED  MODEL  FOR  ANALYSIS 

1 .  General  Description  of  the  Near  Wake 

The  flow  In  the  near  wake  of  a  slender  wedge  is  shown  in  figure  (1). 

At  the  vertex  the  unlfo*-  supersonic  stream  on  each  side  of  the  wedge  is  turn¬ 
ed  through  an  angle  d  ,  through  an  oblique  shock.  Boundary  layers  are  formed 
on  the  top  and  bottom  surfaces  of  the  wedge.  At  the  base  of  the  wedge,  these 
boundary  layers  expand  rapidly,  and  in  a  complex  manner,  to  a  low  pressure. 

The  separation  point  is  believed  to  occur  on  the  base  at  a  small  distance  from 
2,26 

the  corner.  For  low  supersonic  Mach  numbers,  the  acceleration  and  pres¬ 
sure  gradient  are  very  large  in  a  relatively  small  region  near  the  separation 
point.  Hence  it  is  expected  that  the  expansion  process  may  be  treated  as  an 
isentropic  expansion  and  the  flow  may  be  regarded  approximately  as  an  invis- 
cid  rotational  flow.  In  order  that  the  no -slip  condition  along  the  wedge  surface 
be  satisfied  as  the  flow  approaches  the  trailing  edge,  there  exists  a  viscous 
sublayer,  having  steep  velocity  gradient  and  with  thickness  small  compared  to 

the  boundary -layer  thickness.  The  effect  of  the  sublayer  on  the  outer  part  of 

27 

boundary  layer  is  small  .  Since  the  flow  is  actually  overexpanded  and  then 

2 

recompressed  by  a  lip  shock,  deeply  embedded  within  the  shear  layer  ,  the 
expansion  processes  at  the  trailing  edge  are  further  complicated.  The  strength 
of  the  lip  shock  is  small  only  if  the  ratio  of  the  base  pressure  to  wedge  surface 
pressure  is  not  too  small.  Therefore,  for  flow  over  a  slender  wedge  at  low 
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supersonic  Mach  numbers  the  lip  shock  is  rather  weak  and  the  accompanying 

3 

entropy  changes  may  be  neglected  . 

The  expanded  boundary -layer  flow,  separating  from  the  wedge,  mixes 

under  nearly  constant  pressure  with  the  relatively  low -speed  recirculating 
2  4 

flow  ’  .  Viscous  effects  are  important  primarily  near  the  dividing  stream - 
line,  where  the  velocity  gradient  is  large.  For  most  of  the  distance  from  the 
base  to  the  rear  stagnation  point,  the  dividing  streamline  has  small  curvature. 
The  shear  layer  remains  rather  thin  provided,  again,  that  the  wedge  is  slender 
and  die  Mach  number  is  not  very  large.  Eventually  the  shear -layer  flow  ap¬ 
proaches  the  wake  center  line  at  the  rear  jtagnation  point,  which  is  located  at 

2  \ 

a  distance  of  a  few  body  heights  downstream  of  the  base  of  the  wedge  ’  . 

In  the  vicinity  of  the  rear  stagnation  point,  the  flow  outside  die  divid¬ 
ing  streamline  is  turned  to  a  direction  approximately  parallel  to  the  wake 
center  line.  This  fluid  has  momentum  larger  than  the  momentum  at  the  divid¬ 
ing  streamline,  and  therefore  moves  downstream  to  form  the  far  wake.  A 
compression  accompanies  the  turning  of  the  flow,  and  the  compression  waves 
coalesce  to  form  the  wake  shock.  A  fluid  particle  along  the  dividing  stream¬ 
line  is  brought  to  rest  at  the  rear  stagnation  point,  since  it  has  a  total  pres  - 
sure  just  large  enough  to  sustain  the  recompression.  The  fluid  with  lower 
momentum  turns  back  to  the  base  region  and  recirculates.  The  velocity  of  the 
recirculating  flow  is  relatively  small  compared  with  the  velocity  of  the  external 
flow.  The  shear  force  acting  or  the  recirculating  flow  along  the  dividing  stream- 
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line  produces  an  accelerating  torque,  whereas  the  shear  force  acting  along  the 
base  of  the  wedge  produces  a  retarding  torque.  The  net  effect  of  the  acceler¬ 
ating  and  retarding  torques  is  to  allow  a  steady  motion  in  the  recirculation 
region. 

2 .  Motivation,  Assumptions  and  Description  of  the  Present  Model . 

As  suggested  before,  the  near -wake  flow  comprises  three  flow  regions 
of  different  physical  character,  each  described  by  approximate  flow  equations 
of  a  different  mathematical  type.  These  regions  are  the  inviscid  external  flow, 
the  viscous -layer  flow,  and  the  recirculating  flow,  and  are  described  approxi¬ 
mately  by  pure  hyperbolic  equations,  parabolic  equations  and  elliptic  equations 
.respectively.  A  description  of  the  interactions  between  the  regions  will  be  an 
'  essential  part  of  any  analysis  of  the  near -wake  flow  field.  However,  it  seems 
clear  that  the  solution  of  the  problem  depends  largely  on  the  success  in  formu¬ 
lating  essentially  separate  problems  for  the  different  regions  which  may  be 
studied  independently,  and  then  using  some  matching  procedure  to  obtain  a 
unique  solution  for  the  whole  near -wake  flow. 

23 

In  the  study  of  the  recirculating  flow  the  idea  proposed  by  Batchelor 

24 

for  incompressible  laminar  wakes  is  adopted,  in  the  manner  of  Lenard  and 

25 

Guderley  and  Greene  .  That  is,  an  attempt  is  made  to  obtain  an  asymptotic 
solution  for  the  recirculation  region  which  is  valid  in  the  limit  as  the  Reynolds 
number  based  on  external  flow  properties  tends  to  infinity,  with  the  external 
Mach  number  held  fixed.  Because  the  Reynolds  number  is  large  the  viscous 
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forces  acting  on  the  fluid  are  small  everywhere  except  in  the  neighborhood  of 
certain  surfaces  where  the  velocity  gradient  is  large  These  surfaces  enclose 
the  recirculating  flow.  The  region  where  the  viscous  effects  are  not  negligible 
becomes  smaller  and  smaller  as  the  Reynolds  number  increases.  Thus,  in  the 
limit,  the  flow  in  the  interior  of  the  recirculation  region  may  be  regarded  as 
inviscid.  The  solution,  of  course,  is  to  be  used  not  at  arbitrarily  large  Reynolds 

numbers  but  at  Reynolds  numbers  such  that  the  flow  actually  is  laminar . 

2,  4 

From  experimental  evidence  ’  ,  it  is  clear  that  the  length  of  the  vis¬ 
cous  layer  which  accelerates  the  fluid  in  the  inviscid  core  is  greater  than  the 

length  of  the  viscous  layer  which  decelerates  this  fluid.  It  is  expected,  as  pro- 
23 

posed  by  Batchelor  ,  that  as  the  free -stream  Reynolds  number  becomes  large 
the  velocity  in  the  recirculation  region  is  of  the  same  order  as  the  velocity  in 
the  external  flow,  although  numerically  smaller.  In  the  limit,  therefore,  the 
Reynolds  number  based  on  the  typical  properties  for  the  recirculation  region 
will  also  approach  infinity,  although  in  any  application  the  numerical  value  will 
be  considerably  smaller  than  the  Reynolds  number  based  on  external  flow  prop¬ 
erties  .  However,  it  is  believed  that  the  results  will  provide  a  useful  approxi- 
me'ion  at  the  larger  values  of  Reynolds  number  for  which  the  flow  remains 
laminar. 

Since  experiments  suggest  that  the  thermal  energy  of  the  fluid  in  the 
recirculating  flow  is  of  the  same  order  as  the  kinetic  energy  in  the  external 
flow,  which  is  numerically  considerably  larger  than  the  kinetic  energy  of  the 
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recirculating  flow,  the  Mach  number  based  on  typical  properties  for  the  recir¬ 
culating  flow  may  be  assumed  small.  This  assumption  implies  that  the  recir¬ 
culation  region  may  be  treated  approximately  as  an  incompressible  rotational 
flow.  The  effects  of  nonzero  Mach  number  may  be  shown  to  be  rather  small. 

2( 

If  the  energy  and  momentum  equations  are  operated  on  in  the  manner  of  Prandtl 
22 

and  Batchelor  ,  one  may  obtain  the  result  that  the  flow  in  the  interior  of  the 
recirculation  region  is  inviscid  and  incompressible  with  constant  vorticity  and 
temperature . 

In  order  to  study  the  complex  flow  outside  the  dividing  streamline  cer¬ 
tain  simplifications  are  needed  to  provide  a  simple  but  meaningful  solution. 

The  most  important  omissions  are  the  details  of  the  flow  near  the  trailing-edge 
stagnation  point,  where  the  lip  shock  is  formed,  and  the  details  of  the  flow  near 
the  rear  stagnation  point,  where  the  wake  shock  originates.  Experiments  sug¬ 
gest  that  for  a  slender  wedge  at  low  supersonic  Mach  number  the  lip  shock  is 
3 

rather  weak  ,  and  in  the  limit  of  large  Reynolds  number  the  dimensions  of 
these  omitted  regions  are  of  higher  order  than  the  wake  length.  Therefore, (the 
expansion  near  the  trailing-edge  stagnation  point  and  the  recompression  near 
the  rear  stagnation  point  may  be  assumed  to  be  shock-free  and  may  be  replaced 
by  isentropic  processes  .  It  is  also  assumed  that  the  dividing  streamline  is 
straight  and  that  the  pressure  is  nearly  constant  except  in  the  omitted  regions . 

Finally,  for  simplicity  it  is  assumed  that  there  is  no  heat  transfer  to 
the  wedge  surfaces;  the  Prandtl  number  is  unity;  the  product  of  density  and 
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viscosity  is  constant  throughout  the  viscous  layers;  and  the  velocity  profile 
prior  to  the  corner  expansion  is  a  Blasius  velocity  profile  in  suitable  trans¬ 
formed  coordinates. 

Based  on  the  simplifications  and  approximations  mentioned  above,  the 
configuration  of  the  present  model  of  the  near-wake  flow  field  is  shown  in  fig¬ 
ure  (2).  The  thin  boundary -layer  flow  separating  from  the  trailing  edge  of  the 
wedge  expands  rapidly  and  isentropically  to  a  base  pressure .  This  expanded 
boundary -layer  flow  then  mixes  under  nearly  constant  pressure  with  the  rela¬ 
tively  low-speed  recirculating  flow,  confined  in  a  right  triangular  region  with  a 
wake  angle  0Q  between  the  dividing  streamline  and  the  wake  center  line.  The 
mixing  is  limited  to  a  narrow  band  around  the  dividing  streamline  which  has  a 
negligible  curvature.  Eventually  this  shear-layer  flow  reaches  the  wake  center 
line,  and  as  it  begins  to  turn  compression  waves  are  emitted.  The  low-momen¬ 
tum  fluid  below  the  dividing  streamline  is  turned  back  and  recirculates  in  the 
recirculation  region.  The  high-momentum  fluid  above  the  dividing  streamline 
escapes  downstream  and  constitutes  the  far  wake  region.  The  compression  con¬ 
tinues  until  the  external  flow  turns  parallel  to  the  wake  center  line .  The  flow 
in  the  right  triangular  region  may  be  regarded  as  incompressible  and  inviscid 
with  constant  vorticity  and  temperature .  Since  the  no-slip  condition  must  be 
satisfied  at  the  base,  a  base  boundary  layer  is  formed  at  the  base  of  the  wedge. 
This  base  boundary  layer  is  thin  in  the  limit  of  large  Reynolds  number,  and 
produces  resistance  to  the  recirculating  flow  against  the  driving  force  of  the 
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shear  layer. 

3.  Criterion  for  a  Unique  Near  Wake  Solution 

Since  a  fluid  particle  in  the  recirculation  region  undergoes  a  laminar 
steady  motion  along  closed  streamlines,  the  net  gain  in  angular  momentum  of 
each  fluid  particle,  as  it  moves  through  a  complete  cycle,  must  be  zero.  This 
condition  should  be  true  for  all  particles  in  either  the  inviscid  core  flow  or  in 
the  continuous  viscous  layers  enclosing  the  recirculating  flow.  When  a  stream¬ 
line  is  taken  as  a  control  surface,  then  the  abcve  condition  implies  that  the 
angular  momentum  of  the  fluid  bounded  by  this  control  surface  must  be  con¬ 
served.  For  convenience  the  control  surface  may  be  t  v,osen  as  the  streamline 
enclosing  the  whole  recirculation  region,  which  consists  of  the  dividing  stream¬ 
line,  the  wake  center  line  and  the  base  of  the  wedge .  The  shear  along  the  di¬ 
viding  streamline  produces  an  accelerating  torque  whereas  the  shear  along  the 
base  of  wedge  produces  a  retarding  torque.  Since  the  wake  behind  a  symmetric 
wedge  is  symmetric  with  respect  to  the  wake  center  line,  there  is  no  shear 
along  the  wake  center  line.  Therefore,  the  accelerating  torque  from  the  shear 
layer  should  always  be  equal  to  the  retarding  torque  from  the  base  boundary 
layer,  because  there  should  be  no  change  in  total  angular  momentum  of  the  re¬ 
circulating  flow.  This  implies  that  for  a  given  flow  condition  before  the  trail¬ 
ing  edge  expansion  and  for  a  given  wake  angle  there  exists  an  equilibrium  vortidty 
in  the  core  flow  so  as  to  satisfy  the  requirement  of  conservation  of  angular  mo- 
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mentum  for  the  recirculation  region . 

When  the  shear-layer  flow  approaches  the  wake  center  line  it  starts 

to  turn  and  the  pressure  begins  to  increase.  The  recompression  has  been 

assumed  to  take  place  in  a  region  having  dimensions  of  higher  order  than  the 

wake  length.  In  this  region  the  pressure  and  inertia  effects  are  much  larger 

than  the  viscous  effects,  and  so  the  flow  near  and  on  the  dividing  streamline 

may  be  assumed  to  be  recompressed  approximately  isentropically .  A  fluid 

particle  along  the  dividing  streamline  is  brought  isentropically  from  a  large 

velocity,  corresponding  to  the  condition  prior  to  the  recompression,  to  rest 

at  the  rear  stagnation  point.  Accordingly,  its  static  pressure  increases  from 

the  base  pressure  to  the  stagnation  pressure .  Experimental  results  show  that 

2 

at  the  rear  stagnation  pointthe  flow  is  partially  recompres9ed  ,and  he  turning  of 
the  external  flow  to  the  direction  parallel  to  the  wake  center  line  is  not  yet  com¬ 
pleted.  Hence,  further  recompression  is  required  downstream  of  the  rear 
stagnation  point  until  the  external  flow  is  parallel  to  the  wake  center  line. 

If  a  uniform  supersonic  flow  is  turned  through  a  given  angle  from  a 
specified  initial  static  pressure,  the  final  value  of  the  pressure  depends  on  the 

initial  value  of  the  Mach  number,  and  it  is  known  that  the  final  pressure  has  a 

28 

minimum  value  for  a  particular  initial  Mach  number.  Near  the  rear  stagnation 
point,  for  a  streamline  just  outside  the  dividing  streamline  ,  the  velocity  of  a 
fluid  particle  is  decreased  and  the  flow  is  turned  to  a  direction  almost  parallel 
to  the  wake  center  line .  Since  the  turning  actually  is  completed  downstream  of 
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the  rear  stagnation  point,  the  pressure  rise  would  be  overestimated  if  it  were 

assumed  that  the  streamline  becomes  parallel  to  the  center  line  at  a  point  very 

close  to  the  rear  stagnation  point.  To  compensate  partially  for  this  error,  the 

streamline  having  a  minimum  pressure  rise  is  considered.  It  is  assumed  that 

this  streamline  is  turned  through  a  simple  wave  to  a  direction  parallel  to  the 

wake  center  line,  and  that  the  pressure  thus  attained  is  equal  to  the  pressure 

at  the  rear  stagnation  point.  This  requirement  will  be  referred  to  as  the  re- 

9 

attachment  condition.  It  is  similar  to  Chapman's  criterion  of  reattachment 
but  is  based  on  flow  conditions  for  a  particular  streamline  inside  the  shear 
layer  rather  than  outside  the  entire  shear  layer.  Since  for  a  given  flow  con¬ 
dition  before  the  trailing-edge  expansion,  the  shear -layer  velocity  distribution 
is  dependent  both  on  the  vorticity  of  the  inviscid  core  flow  and  on  the  wake  angle, 
the  reattachment  condition  must  also  be  dependent  on  the  same  factors . 

For  a  given  flow  condition  before  the  trailing-edge  expansion,  the  shear - 
layer  velocity  distribution  varies  with  the  change  in  both  the  inviscid  flow  velo¬ 
city  at  the  inner  edge  of  the  shear  layer  and  the  external -flow  velocity;  in  other 
words,  it  depends  on  the  vorticity  of  the  inviscid  core  flow  and  on  the  wake 
angle .  The  velocity  distribution  in  the  base  boundary  layer  is  determined  by 
the  velocity  distribution  at  the  outer  edge  of  the  base  boundary  layer  which  is, 
in  turn,  a  function  of  the  wake  angle  and  the  vorticity  of  the  inviscid  core  flow. 
Hence,  if  the  wake  angle  and  the  vorticity  of  the  inviscid  core  flow  are  obtain¬ 
able,  the  near-wake  flow  properties  may  be  determined.  Therefore,  two  inde- 
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pendent  conditions,  which  involve  both  the  wake  angle  and  the  vorticity  of  the 
inviscid  core  flow,  are  needed  to  provide  a  unique  solution  for  the  near-wake 
flow .  These  independent  conditions  are  the  conservation  of  angular  momentum 
for  the  recirculation  region,  and  die  reattachment  condition  at  the  rear  stagna¬ 
tion  point .  It  should  be  pointed  out  that  the  reattachment  condition  is  considered 
ar  tentative,  since  it  is  not  derived  from  the  solutions  obtained.  In  the  present 
context,  solutions  for  the  flow  details  near  separation  and  near  the  rear  stagna¬ 
tion  point  are  regarded  as  part  of  a  higher  approximation.  A  solution  valid  near 
the  rear  stagnation  point,  and  perhaps  other  parts  of  the  next  approximation, 
would  have  to  be  studied  in  order  to  devise  some  sort  of  reattachment  condition 
in  a  systematic  way.  Presumably  it  could  then  also  be  shown  that  the  error 
really  is  of  higher  order  in  the  present  formulation  of  the  condition  that  angular 


momentum  must  be  conserved . 


Ill 


GENERAL  RESULTS  FOR  RECIRCULATING  FLOW  WITH  LARGE 
LOCAL  REYNOLDS  NUMBER  AND  SMALL  LOCAL  MACH  NUMBER 

1  General  Remarks  and  Formulation 


20  21  22  23 

The  arguments  given  by  Prandtl  ’  and  by  Batchelor  '  for  incom¬ 
pressible  flow  are  extended  here  for  the  steady  flow  at  large  local  Reynolds 

number  in  the  confined  region  behind  a  blunt-based  body  immersed  in  a  uni- 

24 

form  supersonic  stream.  The  derivation  is  also  an  extension  of  Lenard's 
work,  because  an  estimate  is  obtained  to  demonstrate  that  the  effect  of  nonzero 
Mach  number  in  the  recirculating  flow  is  probably  quite  small.  Diffusive  ex¬ 
changes  are  supposed  to  be  so  small  that  they  may  be  disregarded  wherever 
there  areno  large  velocity  and  temperature  gradients  or  accumulative  effects  . 
That  is,  for  steady  laminar  motion  along  closed  streamlines  at  large  local 
Reynolds  number,  viscous  forces  acting  on  a  fluid  element  and  thermal  ex¬ 
changes  between  the  element  and  its  surroundings  are  small  everywhere  except 
in  thin  viscous  layers  where  large  gradients  in  flow  properties  are  encountered. 
Such  thin  layers  enclose  the  recirculating  flow  behind  the  blunt-based  body. 
Hence,  for  large  local  Reynolds  number,  the  flow  in  the  interior  of  the  recircu¬ 
lation  region  may  essentially  be  considered  as  inviscid  and  nonconducting.  The 
assumption  of  small  diffusive  exchanges  implies  that  thenondimensional  forms 
of  the  viscous  forces  and  of  the  thermal  exchanges  are  small  compared  with 
unity,  provided  that  the  characteristic  values  used  in  the  definition  of  local 
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Reynolds  number  and  local  Mach  number  for  the  recirculating  flow  are  chosen 

as  typical  values  for  nondimens ionalization. 

For  the  present  model,  it  is  assumed  that  the  Reynolds  number  Ree* 

based  on  the  external-flow  conditions  and  the  wake  length,  is  large,  and  the 

flow  is  lammar.  Thus  as  Re^  approaches  infinity  as  a  limit,  the  thickness  of  the 

laminar  viscous  layers,  at  the  wedge  surfaces  and  the  boundary  surfaces  of  the 

wake,  is  calculable  and  is  of  order  Re  In  the  limit  Re  "*  °°,  these  viscous 

e  e 

layers  reduce  to  a  number  of  singular  stream  surfaces  of  the  fluid.  Across  such 
singular  stream  surfaces  the  flow  properties  may  be  discontinuous.  For  the  lam- 
inarwakeof  a  wedge,  the  portionof  thewake  boundary  which  yields  resistance  to 
the  recirculating  flow  is  shorter  than  the  portionof  the  boundary  which  produces 
driving  forces,  and  thevorticity  of  the  inviscid  core  flow  must  take  up  an  equi¬ 
librium  value  such  that  the  flow  is  steady.  Therefore,  the  typical  velocity  of 
the  recirculating  flow  must  be  fairly  high.  In  the  limit  of  large  Re^,  the  typi¬ 
cal  velocity  of  the  recirculating  flow  must  be  of  the  same  order  as  that  of  the 
external  flow,  although  numerically  it  is  smaller.  Since  the  flow  properties  in 
the  recirculation  region  differ  from  those  of  the  external  flow  primarily  through 
factors  involving  powers  of  the  Mach  number  M^,  based  on  external  flow  proper¬ 
ties,  then  the  Reynolds  number  Rer,  based  on  the  typical  flow  properties  of  the 
recirculating  flow,  differs  from  Ree  through  factors  involving  powers  of  M£. 

In  the  limit  Reg  *  with  Mg  held  fixed,  Re^  will  also  approach  infinity. 

Let  the  typical  length,  velocity  magnitude,  temperature,  density,  pres- 
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sure  and  total  enthalpy  for  the  recirculating  flow  be  denoted  by  q^,  T  , 
and  =  c^  respectively .  The  coefficients  of  viscosity  and  thermal 
conductivity  corresponding  to  the  typical  temperature  are  p^  and  respec¬ 
tively.  The  following  is  a  list  of  nondimensional  variables  for  the  recirculation 
region: 


where  M  = 


4  yRT 


=  reference  Mach  number  for  the  recirculation  region. 


The  bar  above  any  symbol  denotes  a  dimensional  quantity.  The  first  and  sec¬ 
ond  viscosity  coefficients  are  denoted  by  p  and  p  *  respectively,  and  S  is  the 
specific  entropy.  The  spatial  coordinates  are  made  nondimensional  by  the  typi¬ 
cal  length  which  is  taken  equal  to  the  base  height.  Later  it  will  be  clear  that 
the  obvious  choice  for  q^  is  |  to^L^/2  |  ,  where  is  the  reference  vorticity 
in  the  inviscid  core  flow.  Also  T^  is  the  stagnation  temperature  of  the  free 

stream,  p  is  the  pressure  after  the  trailing  edge  expansion,  and  P  is  the  den- 

_  c  fT 

sity  corresponding  to  Tf.  It  is  assumed  that  the  Prandtl  number  Pr  =  -r  is 

equal  to  one  and  the  specific  heats  are  constants.  In  a  later  part  of  the  deriva¬ 
tion  the  coefficients  of  viscosity  and  conductivity  will  be  taken  to  be  linear  func¬ 
tions  of  temperature,  and  so  p^  and  k^  are  proportional  to  T  . 


The  nondimensional  equations  which  describe  the  laminar  motion  can 
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be  derived  fron  the  general  form  of  those  equations  given  by  Lagerstrom  : 


7  (pq)  =  0 


(IIM) 


P  q  X  (V  X  q  )  = 


~~~2  7P  +  P7  ^  |7(pVq  ) 

YM  r 


-  „2 


+  27  [p(7  •  q  )  ]  +  7(q  •  7p)  -  q  7  M 


+  (7M)X  (7Xq)  -(7  •  q  )(7p)  -  7  X  7  X  (jiJ  )  ]- 


(III  *2) 


-  If  (Y“I)M  r  r  _  _ 

Pq  VH  =  Ri~  [V  <M7H)J+  — R-—  '  q(7  q) 


+  M7(iq2)+P(7  X  J)X  q  ]}  (III-3) 

y  r  i  I  y(y-DmJ  /  _  2 

pq  VS  =  JL  — —  {$<Vq> 


+  ^.72(q-q)+^7-^(7Xq)xq 

-  2q  •  7(7  •  q  )+  |  7  X  q  |2  j  |  (III-4) 


£L  =  *s 


(111*5) 


i 
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PrW 

where  Re  =  _ -  =  Reynolds  number  for  the  recirculation  region, 

r  Mr  S 

r 

AS  =  S  -  S  ,  and  S  =  -  =  the  typical  entropy  for  the  recirculation  region 

-r  -  Cy  -  4  2 

corresponding  to  Tf,  and  .  When  Reg  is  0(10  ),  then  Re^  is  0(10  ),  as 

will  be  shown  later  when  the  numerical  results  are  presented. 

20  22 

In  the  limit  Re^  -•  ®  ,  following  Prandtl  and  Batchelor  ,  the  non- 
dimensional  equations  (III-2)  to  (III -4)  valid  for  the  interior  of  the  region  of 
recirculating  flow  become 

1  2 

pq  x  u  q  =  — 9  vp  +  pv  %-  (in- 6) 

yM 


pq  •  VH  =  0  ( III  -  7) 

pq  -VS  =  0  (III-8) 

Thus  in  the  interior  of  the  recirculation  region  the  flow  is  regarded  as  inviscid 

and  nonconducting.  However,  since  the  gradients  in  flow  properties  are  large 

near  the  boundary  surfaces  which  enclose  the  recirculating  flow,  the  viscous 

effects  are  important  in  thin  layers  adjacent  to  these  surfaces .  As  the  local 

Reynolds  number  grows,  these  layers  become  thinner,  and  the  typical  length  in 

the  direction  normal  to  the  flow  direction  is  not  L  but  should  instead  be  of  the 

r 

order  L^/Ke^ .  Therefore  in  the  limit  of  large  local  Reynolds  number  theappro- 
priate  viscous -stress  terms  may  be  retained  and  these  regions  are  described 
by  boundary -layer  equations . 

Nondimens ional  orthogonal  streamline  coordinates  £  and  0  are  now 
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’ntroduced,  where  £  is  measured  in  the  streamline  direction  and  ip  is  measured 
normal  to  the  streamline  direction.  The  elements  of  length  in  the  direction  of 
increasing  £,  ip  areh^d^,  h^dj/j  respectively. 

The  integral  conditions  which  arise  from  the  effect  of  diffusive  ex¬ 
changes,  and  which  are  valid  no  matter  how  small  the  coefficients  of  viscosity 
and  thermal  conductivity  may  be,  can  be  derived  from  the  energy  equation  and 
the  momentum  equation.  The  operations  applied  to  these  equations  should  be 
performed  in  such  a  way  that  the  contributions  from  all  terms  other  than  the 
terms  involving  diffusive  exchanges  vanish  identically.  Such  an  operation  is  to 
take  the  line  integral  around  a  closed  streamline  in  the  two-dimensional  flow. 

Since  experimental  results  show  that  M  is  fairly  small,  it  will  be 

2 

assumed  that  the  flow  variables  can  be  expanded  in  power  series  of  M^,  as  in 

30 

the  Janzen -Rayleigh  method  (e.g.,  van  Dyke  ,  1964): 


q 

T 

P 


-  .-2  - 
~  q  +.  M  q  ,  +  ... 
o  r  ^  1 

~  T  +  M2T,  +  ... 
o  r  1 

~P0+  MrPl+  ••• 


2 

p  ~  p0  +  Mr  P1  +  ' 

H~  H  +  M2H,  +  ... 
o  r  1 

2 

p  ~  p  +  M  p,  +  ... 
o  r  1 

2 

pr  ~  pr  +  M  p,  +  ... 
ro  r  rl 

The  effects  of  compressibility  may,  in  principle,  be  studied  by  perturbing  a 
basic  solution  for  incompressible  flow.  The  first  order  corrections  are  pro- 


portional  to  the  square  of  the  typical  local  Mach  number  M  ,  and  higher  approx - 

2 

imations  proceed  by  successive  powers  of  . 

2.  Integral  Condition  Obtained  from  Energy  Equation. 

The  integral  condition  arising  from  the  distribution  of  thermal  energy 
within  the  recirculating  flow  is  now  considered.  The  energy  equation  (III -3) 
may  be  rewritten  in  the  following  form: 

V  7  H  =(pLe^pq)  {v  <M™>  +  (y-lJM^  V  •  [  p'  q  (  V  •  q  )  +  p(q  •  V)  q 

(HI-9) 

where  i^  is  a  unit  vector  along  a  streamline  and  q  may  be  replaced  by  qi^. 
Integration  of  equation  (III -9)  around  a  closed  streamline  inside  the  recircu¬ 
lation  region  gives 

$  [  V  7H  >id«  =(r^(5t){7  dVH  +  <v-')Mr2  7  '  [>'?  •?) 

M(?  •*>?]}  hjd*  (III -10) 

The  left-hand  side  of  this  equation  may  be  written  as<j)|^VH  j  ■  hjd£  ,  which 
vanishes  because  H  is  a  single -valued  function  of  position.  Then  the  exact  in¬ 
tegral  condition  to  be  satisfied  for  every  closed  streamline  is 

^q)[7'<,,7H)]hld«  ’  -<Y-l)Mr2^){7  [M'J  <7  ?) 
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Substitution  of  the  expansion  in  powers  of  into  this  equation  gives 


w]v« 


+  [Mr2]^ 


^0%  ' 


VHi)  +  *Vho>-(^  +0^)’-<Vho> 


tijd?  =  ckm’) 


+  (Y'l)V  LM'q0<''  V7)q0 

J 

(III -12) 

Therefore  the  integral  condition  from  the  energy  equation  may  be  written  as 


* 

j (p-^)  ['■  VH0>  Jhld«  ■  °<Mr)  <m-13> 

when  Ref  is  large,  the  energy  equation^  III -7)  for  the  inviscid  core  flow  can  be 
reduced  to 


H  =  H(  j/j  )  =  H0<4>)+  Mr2H1(<<))+  ... 


(Ill- 14) 


This  means  that  total  enthalpy  is  constant  along  a  streamline.  If  equation  (III-14) 
is  substituted  into  equation  (III— 13),  the  resulting  integral  condition  from  the 
equation  for  the  inviscid  core  may  be  rewritten  as 


f(po"oXhiO 


_5  /  hA 

\MCf  h2  d0 


h:d$  =  CKM^  (III-15) 


From  the  definition  of  the  stream  function  it  follows  that 
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°0q0h2 


1 


Hence  the  operator 
becomes 


A 

dip 


may  be  taken  outside  the  integral,  and  equation  (III -15) 


_d_ 

di/j 


0oVoV5 


0 


Integration  with  respect  to  ip  gives 

D0q0M0hld^  =  constant  (III-16) 

Equation  (III -16)  should  be  valid  for  any  streamline  lying  entirely  in  the  invis- 
cid  core  flow,  including  streamlines  arbitrarily  close  to  the  center  of  rotation, 
where  qQ  =  0.  This  indicates  that  the  constant  in  equation  (III  16)  can  only  be 
zero .  Because 


dHo<0) 

dip 


$ 


f°oVohid«  *  0 

except  for  q^  =  0,  equation  (III -16)  becomes 

dH  (iM 

=  0 
u  Ip 


(III - 1 7) 


That  is,  HQ  = 


*  A 


constant,  or 
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Tq  =  constant  (III- 18) 

From  the  definition  of  T  ,  one  gets  Trt  =  1  •  This  condition  states  that  in  the 

limit  Re^  -  *  the  contribution  from  the  diffusive  exchanges  to  the  rate  of 

change  of  enthalpy  in  the  inviscid  core  of  the  recirculating  flow  must  be  zero 

2 

when  terms  of  order  are  neglected. 

The  existence  of  a  constant  temperature  region  in  the  wake  of  a  body 

is  evident  both  from  theoretical  analysis  and  from  experimental  results . 

31 

Burggraf  ,  in  a  study  of  the  structure  of  the  two-dimensional  flow  inside  a 
square,  has  solved  the  Navier -Stokes  equations  numerically  for  a  constant - 
density  fluid.  The  highest  Reynolds  number  for  which  results  are  presented 
is  400,  which  is  comparable  to  values  anticipated  for  Ref  in  the  present  analy¬ 
sis  Both  vorticity  and  temperature  are  found  to  be  approximately  constant  in 
an  inviscid  core  flow.  The  result  of  constant  temperature  in  the  inviscid  core 
is  obtained  by  taking  the  integral  of  the  constant-density  energy  equation  around 
a  closed  streamline.  Todisco  and  Pallone^  have  made  near -wake  flow  field 
measurements,  and  found  the  existence  of  a  temperature  plateau  in  the  recircu¬ 
lating  flow . 

3.  Integral  Condition  Obtained  from  Momentum  Equation 

The  integral  condition  arising  from  the  velocity  distribution  may  be 
obtained  by  integrating  the  momentum  equation  around  a  closed  streamline  in¬ 
side  the  recirculating  flow.  For  convenience  the  momentum  equation  (III— 2) 
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there  results: 
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(rt)( $r0  ['  x  ^0  < v  x  «o> +  ■  V  +  < 7  •  *o  ’)?0]  -h, d« 

+  (Mr^[  7  X  Ml(7  X  q>  7(M1  7'  <><7-VN 

(P  ,\  -  /P,\  ,  -  /M 

\J0)  7  x  po<7  x  *«o>  \FJ  7V-«o>  \P^j(7'2jlo7)qo 

+  V  X  Mq(V  X  q{)+  ?(M'  V- J1)  +  (7.2m0V)^1:J  -l^d* 


+  0(M  ) 
r 


•feRv!' 


(III-21) 


For  the  interior  of  the  recirculation  region,  an  estimate  of  the  variation  in  en¬ 
tropy  along  a  closed  streamline  can  be  derived  from  equation  (III -4)  and  the 
condition  of  constant  temperature  in  the  Lnviscid  core: 


(III -22) 


Now  the  estimate  of  the  right-hand  side  of  equation  (III -20)  may  be  obtained  by 
substituting  the  estimate  of  the  entropy  variation  from  equation  (III-22)  into  the 
equation  of  state  (III -5) .  Because  equation  (III-5)  can  be  expanded  as 


1_ 

P 


1  +  -  AS  +  . 
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(III -23) 


the  right-hand  side  of  equation  (III -20)  has  the  following  expansion: 


« 


-v*  =°(t^)  v* 


+  +  O  [ (AS)2  ]  (III -24) 


The  first  term  on  the  right-hand  side  of  the  above  equation  is  zero.  By  sub- 

2 

stituting  the  expansions  in  powers  of  and  the  estimate  of  equation  (III -22) 
into  equation  (III -24),  one  finds 


f( 


?)'V«  =  o 


<u‘"fK"‘  ir  ) 


(III- 25) 


Since  T^  =  1,  and  S  s  constant  along  a  streamline,  it  follows  from  die  equa¬ 
tion  of  state  and  the  relation  between  S  and  that  p_  and  p_  are  constant  along 

v  0  0 

P 

streamline.  From  the  component  of  the  momentum  equation  normal  to  stream - 

2 

lines  it  is  found  that  the  changes  in  pressure  normal  to  streamlines  are  CKM^). 
Therefore  p^  and  p^  are  constant  throughout  the  core  region.  Hence  one  ob¬ 


tains 
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Since  (J^  =  Mq(Tq),  =  Mq(Tq),  and  =  constant,  the  resulting  integral 
condition  from  the  momentum  equation  can  be  rewritten,  to  order  one,  as 

$|><w]'hid« 

(III -27) 

The  term  on  the  right-hand  side  vanishes  because  V-q^  ^  0.  Therefore 
$[*  X  (V  X  q0)]-hld€  =0  (III -28) 


(III -26) 


It  is  clear  that  to  order  one  the  recirculating  flow  may  be  considered  as  in¬ 
compressible  because  both  and  p^  are  constants.  The  result  of  equation 
(III- 28)  may,  of  cours<  be  obtained  by  operating  on  the  incompressible  flow 
momentum  equation.  In  order  to  obtain  a  condition  which  is  valid  for  the  en¬ 
tire  inviscid  core,  rather  than  for  a  single  streamline,  equation  (III -6)  is  used. 
By  taking  the  curl,  and  using  V  •  (  7  X  q  )  =  o  and  also  p  =  p(p),  equation  (III -6) 
can  be  reduced  to  the  vorticity  equation  for  inviscid  flow 


(q 


q 


(III-29) 


tO 

r®r 


For  the  two-dimensional  recir- 


where  to  =  7  x  q  and  to  = 
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culating  flow,  the  right-hand  side  is  identically  zero,  because  co  is  always 
— * 

perpendicular  to  q  .  Thus  it  follows  that 

t  •  °<*> 

2 

Let  the  expansion  of  vorticity  in  powers  of  be  defined  by 

2 

CO  ~  Oh  +  M  CO,  +  . .  . 

0  r  1 

and 

0Q(i/))  +  Mr2  0^1/0  +  •  •  • 

hen 

co~p0o0<iM  + (Mr2)  [p1n0(^)+  Po^Oi))]  +  ••• 

hen  the  above  result  is  applied  to  equation  (III -28)  there  results 

X  p0O0<^)"co]'hld^  =°  (III-30) 

— *  — ♦ 

where  i  is  the  unit  vector  normal  to  the  plane  of  the  flow  and  <o  =  co  i 
CO 

Since  Pq  is  a  constant,  equation  (III-30)  can  be  rewritten  as 

£b?no(*>khid«  - 0  (m-31) 

Since  <£>  q^h  d£  i  0  except  when  qQ  =  0,  the  condition  obtained  from  the  mo- 
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mentum  equation  for  the  inviscid  core  flow  can  be  written,  to  order  one,  as 

=  constant  or  =  constant  (111-32) 

From  the  definition  of  |  |  ,  must  be  equal  to  -1.  This  condition  states 

that  in  the  limit  Re  -* »  the  contribution  from  viscous  forces  to  the  rate  of 
r 

change  of  vorticity  in  the  two-dimensional  recirculating  flow  must  be  zero 

2  20  22 
when  terms  of  order  are  neglected .  Prandtl  and  Batchelor  have  found 

that  for  incompressible  flow  with  closed  streamlines  the  vorticity  is  uniformly 
distributed  in  the  inviscid  core  flow.  Since  the  conditions  =  constant  and 
Tg  =  constant  for  the  inviscid  core  flow  are  obtained  by  studying  the  dissipa¬ 
tive  terms,  they  demonstrate  the  fact  that  even  though  the  core  flow  is  inviscid 
and  nonconducting,  and  the  primary  effects  of  diffusive  exchanges  are  limited 
to  the  continuous  viscous  layer  enclosing  the  recirculating  flow,  the  accumula¬ 
tive  effects  of  very  small  diffusive  exchanges  are  experienced  throughout  the 
entire  recirculating  region.  In  other  words,  the  reason  for  the  existence  of 
these  integral  conditions  are  that,  in  the  exact  steady  motion,  the  net  effects  of 
diffusive  exchanges  must  be  exactly  zero.  Therefore,  for  two-dimensional 
flows  with  closed  streamlines,  the  fluid  motion  cannot  be  exactly  steady  until 
the  effects  of  the  small  but  persistent  diffusive  exchanges  have  evened  out  any 
variation  of  temperature  and  vorticity  that  may  have  been  present  initially. 

4.  Discussion  of  Errors  in  the  Incompressible-flow  Approximation  for  the 
Inviscid  Core . 

The  solutions  obtained  in  the  preceding  sections  are  based  on  the 
assumption  that  the  density  of  the  recirculating  fiow  is  approximately  constant. 


-  38  - 


The  errors  resulting  from  this  approximation  depend  on  the  magnitude  of  the 
reference  Mach  number  for  the  recirculation  region,  which  is  known  to  be 
numerically  small.  Accordingly  the  higher -order  terms  in  M^,  which  have 
been  neglected  in  the  previous  incompressible -flow  analysis,  should  be  exam¬ 
ined  so  that  an  estimate  of  the  compressibility  effect  on  the  recirculating  flow 
can  be  obtained. 

2 

If  terms  of  order  are  retained,  the  next  approximation  to  the  inte¬ 
gral  condition  from  the  energy  equation  can  be  obtained  from  equation  (III-3): 


(III -33) 


Since  =  H^(0)  for  the  inviscid  core  flow,  the  left-hand  side  becomes 


where  p^q^h^  =  1.  The  right-hand  side  of  equation  (III -33)  may  be  rearranged 
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where  cc L  =  -  -r-r —  .  Hence  equation  (III-33)  can  be  written  as 

°  hi  2 


2  **0 
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■  hld€ 


(III -34) 


Integrating  equation  (III-34)  with  respect  to  Ip  , 


|[^(Hi(l|,)  +  1r  %)  poqohide 


■  [<v‘)"o]Ifp^  V«d* 


Changing  the  order  of  integration  on  the  right-hand  side  term  and  using 
P0q0h2  =  1# 

2 

)Wld?  =  (v-DWo  A<*)  -(v-i)y(^-)p0q0h,d« 


(III -35) 


-  40  - 


where  A(0)=  i  6  —  —  h  d£d0  =  1  6h0h,  d £ dl/>  .  Hence  A(0)  is  the  area 
inside  the  closed  streamline .  Also 


Since  Stokes'  theorem  gives 


-  woAO P) 


then  equation  (III-35)  can  be  reduced  to 
dHjOfi) 
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(III-36) 


Integrating  the  above  equation  with  respect  to  ip  and  rearranging, 

n  \Or  '  •  2 

q~  \2 


Ti =  *r  %  + 


4>l  T“ 

'l 


J. 

2 


I  ^O^l^ 


dip 


(III-37) 


In  order  to  obtain  the  exact  value  of  T^,  it  would  be  necessary  to  eval¬ 
uate  the  second  term  on  the  right-hand  side  of  equation  (III -37) .  Since  the  ex¬ 
act  functional  dependence  of  h^  is  not  known,  the  value  of  the  integral  can  not 
be  obtained  analytically.  However  an  estimate  of  the  integral  can  be  performed, 
so  that  the  contribution  of  may  be  estimated.  The  term  in  the  curly  bracket 


typical  of  A(  l/) ) 


therefore  the  integral  of  equation  (III -37)  will  be  estimated  as 


AJ/) 

typical 


of  Mil* ) 


All) 

typical  of  A(  ip ) 
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typical  of  A((i> ) 


It  is  assumed  that 
equation  (III-37)  as 


s  1  Hence  an  estimate  of  may  be  obtained  from 


T1  *  "  (*T  )  1%  ‘  Lq0  J  typical  of  A(V)}  (III-38) 
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From  the  solution  for  the  inviscid  core  flow,  which  will  be  treated  in  (IV),  one 

may  show  tha*-  q^  is  always  less  than  unity.  Actually  I  q^l  max  *  0.85,  and 

so  the  quantity  in  the  curly  bracket  of  equation  (III-38)  must  be  always  smaller 

than  unity,  with  maximum  value  probably  less  than  0.5.  Furthermore,  for 

2 

y  =  1.4,  the  numerical  factor  is  equal  to  0.2.  Therefore,  the  term  T 

is  believed  to  be  very  small  in  comparison  with  Tq,  and  for  all  practical  pur- 
2 

pose  T  may  be  considered  as  zero. 

It  is  known  that  in  the  inviscid  core  flow  p^  constant.  Thus  the 
momentum  equation  can  be  written  as 


(III-39) 


(III -40) 


Solving  these  equations 


P!  = 
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Po'b 
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2  '  ‘  VJ 

(III-41) 


where  <l>  =  —  r  ,  and  l(f>  |  is  the  mass  flow  per  unit  time  between  the 

c  \  r  c  _  _ 

center  of  rotation  and  the  enclosing  boundary  surface.  Since  a  9gmax 

amfc'iHlhe  distance  between  the  center  of  rotation  and  the  location  where 


qQ  =  qQ  mflX  ^  approximately  equal  to  L^/2  ,  the  mean  velocity  between 
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these  two  points  will  be  estimated  as  Therefore  ([>£  will  be  estimated  as 
-  0.5.  By  the  definition  of  the  typical  values  of  density,  temperature  and  pres 
sure,  it  is  clear  that  Pq  =  Pq  =  Tq  =  1,  and  so  the  equation  of  state  may 
be  written  as 


P1  X  P1  +  T1 


(III-42) 


Since  T  a  0, 
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By  using  equation  (III -43)  and  the  expansion  of  0  in  powers  of  ,  there  results 

«  =  n  =  -i 
o  o 


or 
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From  the  energy  equation,  die  equation  of  state,  and  die  results  so  far  obtained, 
it  has  been  shown  that  along  a  streamline  the  variation  ir.  entropy  is  AS  = 

*2  v  / 


/  M  \ 
°U) 


and  therefore 


jhj  d£  =  0^  J  .  In  the  expanded  integral 


condition  (III-21)  all  but  one  of  the  terms  of  order  are  zero,  and  so  the  next 
approximation  to  the  condition  becomes 


$[7XV«]  V*  =  0 


(III -44) 


Substituting  the  expression  derived  for  into  equation  (III-44)  one  obtains 
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which  may  be  rewritten  as 
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After  rearranging  and  then  integrating  with  respect  to  tp , 
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the  integral  of  equation  (III-46)  will  be  estimated  as 


f  dtp  (III-46) 
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Then  equation  (III -46)  will  be  estimated  as 


Q 

1 


& 


Y 


typical  in  A(  Ip ) 


co. 


typical  in  A(  ip )  J 


-<*-*> 


* 

» 

% 


Nearing  the  center  of  rotation,  ip  **  lpQ  and  qQ  “*  0,  and  so 


coi  -  0 


Approaching  the  edge  of  the  inviscid  core  flow,  Ip  0  and  0  *  q^  <  q^  max 
s  0.85.  Since  <l'c  &  -  0.5,  |  C0^  |  might  be  as  large  as  0.5  near  the  stagnation 
points.  Elsewhere  in  the  inviscid  core  flow.it  would  be  expected  that  I  co^  |  is 
considerably  smaller  than  0.5.  Thoughout  most  of  the  region,  therefore 
I  co^  |  is  believed  to  be  quite  small  in  comparison  with  |  co^  |  . 


/ 


IV 

THE  INVISCID  CORE  FLOW 
1 .  General  Description  and  Formulation  of  the  Problem 

The  results  from  previous  sections  have  shown  that  when  the  wake 

Reynolds  number  Re^  is  large  and  the  wake  Mach  number  is  small,  the 

20  22 

model  proposed  by  Prandtl  and  Batchelor  should  correctly  describe  the 
interior  of  the  recirculation  region.  Accordingly  two  conditions,  which  im¬ 
ply  uniform  temperature  and  vorticity,  have  been  deducted  for  this  inviscid 
core  flow. 

The  static  pressure  in  the  near  wake  is  known  to  be  nearly  constant, 
except  in  the  vicinity  of  the  trail  Jng-edge  stagnation  point  and  the  rear  stagna¬ 
tion  point.  However,  the  regions  near  these  points  are  assumed  to  be  small  in 
comparison  with  the  length  of  the  wake,  provided  that  the  Mach  number  is  not 
too  large.  Away  from  the  stagnation  points,  the  curvature  of  the  dividing 
streamline  may  be  regarded  as  small.  Furthermore,  since  the  wake  Reynolds 
number  is  large,  the  regions  near  the  enclosing  surfaces,  where  the  diffusive 
exchanges  are  large,  can  be  considered  small.  Consequently,  die  interior  of 
the  recirculation  region,  where  the  recirculating  flow  is  inviscid  and  incom¬ 
pressible  with  uniform  vorticity,  may  be  approximated  as  a  right  triangle 
bounded  by  the  base  of  the  symmetric  wedge,  the  wake  center  line  and  the  di¬ 
viding  streamline . 

It  is  convenient  to  work  with  a  stream  function  and  to  use  cylindrical 
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coordinates  5  and  r,  with  origin  at  the  rear  stagnation  point  and  £  measured 
clockwise  from  die  center  line  of  the  rec^  culation  region  (see  figure  3).  For 
plane  incompressible  flow  the  stream  function  is  defined  by 


1  -  - 
where  ^  and  -  r  ^  are  the  velocity  components  in  the  0  and  r  directions 

respectively.  Therefore  the  governing  differential  equation  is  the  Poisson 


equation 


V 2  =  to 


(IV-l) 


where  co  =  co^  and  is  a  constant.  This  equation  expresses  the  physical  fact 
that  the  flow  is  incompressible  and  inviscid  with  a  uniform  distribution  in  vor- 
ticity.  The  boundary  conditions  to  be  satisfied  are 


(r  cos0  _\ 

-^sy°>  0  )  =  0  (IV-2) 


where  rQ  is  the  length  of  the  dividing  streamline  and  0Q  is  the  wake  angle . 
These  boundary  conditions  indicate  that  in  the  triangular  region  the  mass  is 


conserved;  that  is,  there  is  no  mass  exchange  with  the  surroundings .  In  order 
to  nondimens  ionalize  and  normalize  the  governing  differential  equation  and 
boundary  conditions,  new  variables  are  introduced  in  the  following  manner: 
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Tlien  the  nondimens ional  form  of  equations  (IV-1)  and  (IV-2)  can  be  written  as 


I  _L_ 

r  ar  r2g 2 


aJjL*  =  i 
ae2 


0*(r,O)  =  0*(r,  1) 


cos  6 


(IV-3) 


0  (IV-4) 


2 .  Method  of  Solution 

In  this  problem  a  particular  integral  of  the  differential  equation  (IV-3) 
can  be  obtained  by  defining 


0*( r,  e ;  eo)  =  Mr,  0 ;  0Q)  +  X(r,  0 ;  0o> 


(IV-5) 


where  >&(r,  0 ;  0q)  is  the  particular  integral  and  satisfies  the  following  differ¬ 
ential  equation  and  boundary  conditions: 


ajp  i  ajL  _j_  a^ 

.  2  r  ar  2rf  *fl2 
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=  1 


(IV- 6) 


Mr,  0;  0)  =  M  r,l:  0)  =  0 
o  o 


(IV-7) 


A  solution  of  equations  (IV-6)  and  (IV-7)  is 
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*(r,  0;  eo)  = 
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0)  ]  +  sin  [  20q0  ]  -  sin 
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(IV-8) 


If  equations  (1V-5)  and  (IV-8)  are  substituted  into  equation  (IV-3)  and  boundary 
conditions  (IV-4),  the  differential  equation  becomes 
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and  die  boundary  conditions  reduce  to 


(IV-9) 


X(r,0;  0^)  =  x  (r,  1;  0q)  =  0 


(IV-  10a,  b) 
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(IV-lOc) 


Although  it  may  be  possible  by  standard  methods  to  obtain  a  solution  which  sati- 

fies  (IV-9)  and  (IV-10)  to  any  desired  degree  of  accuracy,  it  appears  that  any 

such  solution  would  be  very  complicated.  Since  the  wake  angle  ?  is  known  to 

c 

be  rather  small,  a  further  simplification  can  be  made  by  expanding  tp  in  a  (low¬ 
er  series  in  0q.  Equation  (IV-8)  has  the  following  series  expansion  in  0q: 


*<r,  0;  0q)  =  0^1(r,0)+  0^  *3(r,  0)  +  0(0^) 


M&r*C*K%  m<rnm 
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<e  -  e  ) 


l]  -4 

-)  +  ^[(e- 


203  + 


eV_i 

0  >  6  j 


+  O<<0 


(IV-11) 


The  solution  for  x  must  satisfy  equations  (IV -9)  and  (IV- 10).  Since 
the  boundary  condition  given  by  equation  (IV-lOc)  may  be  approximated  by  a 

r 

Fourier  sine  series,  the  method  of  separation  of  variables  can  be  employed.'' 
The  largest  term  of  the  solution  is  of  the  form 


.ff/0o 


sin  it  6 


where  7r/0o»  1  because  0  <  0  «  1.  Therefore  the  solution  decays  fairly 
rapidly  as  (1-r)  increases.  Consequently,  the  contribution  of  x(r,  0;  0q)  to  ip* 
may  be  considered  as  limited  to  the  region  where  (1  -r)  iS  small.  Since 

~  -  exp  [-7r(l-r)/0Q  ]  as  r-*  1,  it  appears  that,  for  0q  -*  0,  x  is  a  function 
of  the  stretched  coordinate 


With  this  new  independent  variable,  the  expansion  of  x(r,  8 ;  0q)  is 


X(r,0;eo) 


e2X,<r*.  0)  ^  03  X9(r*,0)  +  04X,(r*,0)  +  0(0^) 
o  1  o  2  oo  o 


(IV-12) 


Hence  die  expansion  of  is 
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+  . . . 


(IV-13) 

When  equation  (IV-12)  is  substituted  into  equation  (IV-9)  and  (IV-lOa, b, c), 
there  results 


Xir*r*+Xlt3  +  «0C*2r,r,-xlr<  +  )(2ee+2r.!(ia(i:1+  .  0 


100 


X](r*  0)+  0v(r*,O)  +  ...  = 


(IV-14) 

(IV-15a) 


Xj(r*,  1)  +  eoX2(r*,l)  +  . . .  - 


(IV-l5b) 


*i(0's>  +  e,[*2<<w+ 


lr* 

2,2 

8  Alr*r* 


«M»>] 


+  82rv(°ie)  +  (kD_  foo.o-e2) 

°  L  3  s  - — -'x2r,(o,e)]+ 


=  l[e-92]+^[-I9+4e2+ie3.,64j  + 


(IV-15C) 


3.  The  Solution  for  the  Inviscid  Core  Flow 

From  equations  (IV-,4)  and  (.V-  15a,  b,  c),  dte  governing  differential 
equation  and  boundary  conditions  for  Xjfr,  9)  can  be  written  as 


*  1 r*r* 


f  X 


100  = 


0 


(IV- 16) 
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Xj(r*.0)  =  Xj(rM)  =  0  •  (IV-17a,b) 

x^o.e)  =  i(e  -  e2)  (iv-i7c) 


The  above  differential  equation  and  boundary  conditions  describe  inviscid  in¬ 
compressible  flow  in  a  semi -infinite  strip  bounded  by  walls  at  0  =  0  and  0=1; 
the  boundary  condition  at  r*  =  0  is  prescribed  by  equation  (IV- 17c);  the  other 
condition  needed  for  the  solution  of  the  above  problem  can  be  obtained  by  the 
requirement  that  the  solution  must  be  bounded  as  r*  “*  ®  .  Since  the  boundary 
condition  at  r*  =  0  has  the  Fourier  series  expansion 


i<0  -  6)  =  A, 


_4_  Y  _ i 

n3  ft  <2n+ir 


sin(2n+l)7r  0 


(IV-18) 


and  since  the  method  of  separation  of  variables  can  be  applied,  the  solution 
for  Xj(r*.  0)  is 


Xj(r *,  0)  =  A  Y.  - ^-3  exp  [-(2n+l)7r  r*  ] 

ff  U  (2n+l) 


—2 

Therefore  the  solution  to  order  0  can  be  written  as 

o 


sin(2n+l)  ir  0 

(IV- 19) 


*'<r-  e;eo>  1  90 


+  0<$ 


2  ® 
r  (ft2  4  V 

2  <e  •  e)+_3  L 
*  0 


[exp[-(2n+lXl-r)?r/0(3] 
L  (2n+l)3 


sin(2n+l)?r0 


(IV- 20) 
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For  small  values  of  the  wake  angle  0q,  this  solution  is  equivalent  to  the  so¬ 
lution  for  the  flow  in  a  circular  sector  with  angle  0  and  radius  r  .  The  flow 

o  o 

in  the  circular  sector  with  constant  vorticity  would  be  described  approximately 
as  follows: 


r(r,  0;ffo)  =  [-  Ty  (0  -  e2)  0o2  +  Xj(r,  0) 


b  xi  i_ 

3r2  +7  V  *02 


X^r.O)  =  Xj(r,  1  )  =  0 


+  0<£> 


XjU-S)  =  He  -  e2) 


and  Xj(r.  0)  would  be  obtained  by  using  the  method  of  separation  of  variables. 
The  solution  is 


r<<>2'e)+  h 

7T 


•  r  C(2n+1)  ir /0q] 
£  (2n+l)3 


sin(2n+l)7T0 


+  0<?3) 


Since  the  contribution  from  the  series  terms  diminishes  as  (1-r)  grows  and 
since  the  powers  of  the  nondimensional  radius  r(2n+^7r/^>  for  (i-r)  small  may 
be  approximated  by  the  exponential  form  e  ^2n+^^  r)ff/®0(  ^g  a^yg  solution 
may  be  approximated  by  equation  (IV-20)  or  vice  versa. 
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(IV-23b) 
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Since  x^  satisfies  homogeneous  boundary  conditions,  x^  satisfies  the  follow¬ 
ing  differential  equation  and  boundary  conditions: 


l2r*r< 


„(c) 

200 


x^'*,.*  +x':;.  =  0 


(IV- 24) 


X(2CV:0)  =  X(?C)(r*,l)  =  0 


(IV25a,  b) 


x(2c)(o,  0)=  -|(i-02)Y  — 
ir  XT  (2n+l) 


sin(2n+l) 7r  0  (IV-25c) 


The  other  boundary  condition  is  the  condition  that  the  solution  must  be  bounded 
as  r*  -*  00  . 

(c) 

The  boundary  condition  for  x2  (r*,0  )  at  r*  =  0  is  expressed  as  the 
product  of  a  polynomial  and  an  infinite  sine  series .  This  sine  series  is  known 
to  converge  uniformly  in  the  region  0  i  0  £  1,  because  the  coefficients  are 
0(l/n2)  asn-*  ® .  Equations  (IV-24)  and  (IV-25)  can  be  solved  by  the  method 
of  separation  of  variables  provided  that  an  appropriate  Fourier  series  repre¬ 
sentation  of  the  boundary  condition  at  r*  =  0  can  be  obtained.  Since  the 
boundary  condition  is  epxressed  by  a  uniformly  convergent  series  which  con¬ 
verges  fairly  rapidly,  only  a  small  number  of  terms  are  needed  to  provide  a 

(c) 

quite  good  approximation.  If  the  sine  series  appearing  in  x2  (0,  0)  is  approxi¬ 


mated  in  this  manner,  a  number  of  ways  might  be  suggested  to  expand  the 
resulting  expression  in  Fourier  series.  However,  for  a  good  fit  with  a  small 


number  of  terms,  with  term-by-term  differentiation  permitted,  a  Fourier  series 

which  converges  rapidly  is  desirable.  The  procedure  adopted  here  will  provide 

3 

a  series  with  coefficients  which  are  0(1  /n  )asn-*  ®  . 

(c) 

For  the  present  purpose  the  boundary  condition  (0.  6)  may  first 
be  approximated  quite  accurately  by  a  fourth -degree  polynomial  in  0: 

X(2C)(0,  6)  e  £  An0n 

0 

Its  coefficients  are  determined  by  solving  a  set  of  five  algebraic  equations 


The  resulting  approximated  boundary  condition  may  be  written  as 
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Consequently  the  problem  of  equations  (IV-24),(IV-25a,b)  and(IV-26) 
can  be  solved  by  the  method  of  separation  of  variables .  The  solution  is 


L4-r+-Jt-a 

jt  0  [\[(  2n+l)7r]  [<2n+l)ir]/ 


+  - =)exp[-(2n+l)7r  r<]sin(2n+ 1)  tt  0 

[<2n+l)ir]7  J 


r4  [(<*’) 


exp  [-2nirr*]  sin  2n  ir  0 


Thus  the  solution  for  X^r*,  0)  is  expressed  by 


(IV- 27) 


i  *  a\  -  ±  V  \(  r*2  7.4 

X„(r  >  ®)  2  '  I"  2  +  3  3 

1  tt  o  L\  2(2n+l;  <*»+!)» 


36  \  * 

+  - - — -  ]  exp[-(2n+l)  t  rl  sin(2n+l)x  0 

<2n+l)V/ 


ntSh''2””* 


]  sin  2nx0 


(IV-28) 


Because  of  the  presence  of  terms  involving  sin  2n  it  0  ,  the  solution  x^X*,  0)  is  not 
symmetric  with  respect  to  the  line  ?  =  $0q.  The  complete  solution  for  the  invis- 


cid  core  flow  may  be  written  bt 


•  /exp[  -(2n+l)(l-r)ir/0 


(2n+l) 


)sin(2n+ 1)  7r  0 


d  ^  J(- 

if  oL \ 


(l-rr  t  7.4 
2(2n+l)20)2  (2n+l)37r3 


+  ^5  5 y^P  [-<2n+l)(l-r)7r/0o  ]jsin(2n+l)7r  0 

+  72 


(IV-29) 
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represented  by  infinite  series.  For  reasons  already  mentioned  these  infinite  - 
series  representations  are  approximated  by  polynomials  in  0.  If  these  polynom¬ 
ial  approximations  are  to  tx  fairly  accurate  at  every  point,  it  is  necessary  to  en' 
sure  that  the  original  series  be  uniformly  convergent.  Since  Xlr*r*  evaluated  at 
r*  =  0  has  the  following  series  expression 

0» 

Xlr«r*<0’8>  '  *  I  (2n+T)  sln<2n+1>*8 
0 

then  Xlr*r*(0,  0)  =  1  0  <  0  <  1 

but  if  0  is  set  equal  to  zero  or  one,  the  series  gives 

xlr*r*(0l0)  =  xlr*r*((U)  =  0 

Obviously  die  series  expansion  of  Xlr*r*(0,  0 )  does  not  converge  uniformly  at 
0  =  0  and  0  =  1.  For  the  third  or  higher  order  radial  derivatives  of  x  j ,  the 
resulting  series  expressions  evaluated  at  r*  =  0  are  no  longer  convergent  series . 
Therefore  the  approximate  polynomial  representation  of  the  boundary  condition 
can  only  be  applied  to  Xlr,(0,  0)  with  sufficient  accuracy. 


V 


BASE  BOUNDARY  LAYER 

1 .  Formulation 

As  the  wake  Reynolds  number  approaches  infinity,  for  a  first  approxi¬ 
mation  the  recirculation  region  can  be  divided  into  the  inviscid  core  flow  and  the 
viscous -layer  flow.  Since  the  wake  Mach  number  is  small,  the  recirculating 
flow  may  be  regarded  as  incompressible.  The  inviscid  core  flow  has  already 
been  studied  in  (IV).  The  viscous  layers  enclosing  the  inviscid  core  flow  form 
a  thin  continuous  layer  for  large  Reynolds  number.  In  order  to  obtain  an  approxi¬ 
mate  solution,  it  is  necessary  to  break  up  this  continuous  layer  into  separate 
sections  which  may  be  solved  approximately.  These  sections  are  the  shear  lay¬ 
er  which  separates  the  external  flow  from  the  recirculating  flow;  the  base  bound - 

23 

ary  layer  along  the  base  of  the  wedge;  and,  as  mentioned  by  Batchelor  ,  the 
viscous  layer  along  the  wake  center  line  between  the  upper  recirculating  flow  and 
the  lower  recirculating  flow  (which  will  not  be  studied  here). 

In  studying  the  base  boundary  layer  near  the  wake  center  line,  the  velo¬ 
city  distribution  of  the  flow  impinging  on  the  base  may  be  approximated  by  the 
velocity  of  the  inviscid  core  flow,  because  as  the  Reynolds  number  grows  the  vis¬ 
cous  layer  along  the  wake  center  line  becomes  thinner,  and  because  the  velocity 
gradient  across  this  layer  decreases  as  the  fluid  moves  toward  the  base.  In  the 
vicinity  of  the  stagnation  point  at  the  center  of  the  base,  the  velocity  components 


-  60  - 


-  61  - 


parallel  and  perpendicular  to  the  base  for  the  inviscid  core  flow  are  linear  func¬ 
tions  of  the  distance  from  the  wake  center  line  and  the  distance  from  the  base  of 

the  wedge  respectively.  This  type  of  flow  is  therefore  a  two-dimensional  stag- 
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nation -point  flow  for  which  the  solution  is  known  ’  .  The  pressure  in  the  in¬ 

viscid  core  flow  is  given  by  the  Bernoulli's  equation  for  rotational  flow.  For  the 
stagnation -point  viscous -layer  flow  the  component  of  the  momentum  equation  in 
the  direction  parallel  to  the  base  is  independent  of  the  component  normal  to  the 
base,  and  is  the  same  as  the  boundary -layer  equation  for  incompressible  flow. 
Since  the  base  is  flat  and  the  boundary -layer  assumptions  are  satisfied  in  the 
viscous  layer  everywhere  along  the  base  except  very  near  the  stagnation  point  or 
the  point  of  separation  just  below  the  corner,  the  governing  differential  equations 
for  this  layer  are  identical  to  the  flat-plate  boundary -layer  equations  for  incom¬ 
pressible  flow: 


au 

5T 


+ 


(v-i) 


-  d% 


(V-2) 


where  U,  V  are  velocity  components  parallel  and  perpendicular  to  the  base  re¬ 
spectively?  X,  Y  are  the  coordinates  parallel  and  perpendicular  to  the  base,  and 
are  measured  from  the  wake  center  line  and  the  base  of  the  wedge  respectively: 
and  q.  is  the  velocity  in  the  inviscid  core  flow  evaluated  at  the  base.  Because 
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the  no-slip  condition  must  be  satisfied  at  the  base,  and  the  flow  velocity  at  the 
outer  edge  of  the  viscous  layer  approaches  the  inviscid-flow  velocity  evaluated 
at  the  base,  the  boundary  conditions  can  be  written  as 
0  =  V  =  0  at  Y  =  0 


U  =  n 


at  the  outer  edge  of  viscous  layer 

(V-3) 


This  system  of  gcverning  differential  equations  for  a  plane,  incom¬ 
pressible,  laminar  boundary  layer  may  be  transformed  by  introducing  the  appro 

35 

priate  form  of  the  Dorodnitsyn  variables  ’  ; 


_ 

Y  =  q  Y 

b 


u  -  i 


% 


W  =  (  +  ^  Y  S-  )/rT 

\%  %  W  %  )  r 


One  then  obtains 


MJ  aw 
ax  +  aY 


0 


(V-4) 


u^  +  w^3(.-u2)  +  ^ 

JX  JY  qb  5Y2 


(V-5) 


where 


d% 

d X 

U 

U 


The  boundary  conditions  are  transformed  to 
W  =  0  at  Y  =  0 

1  as  Y  -*  ® 
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Hence  in  the  Dorodnitsyn  variables  the  boundary  conditions  are  normalized. 

35 

The  method  of  integral  relations  is  employed  to  solve  the  system  of 
transformed  boundary -layer  equations .  This  method  converts  the  system  of 
partial  differential  equations  into  a  s)  3tem  of  ordinary  differential  equations  and 
permits  the  use  of  well -developed  numerical  methods  of  solution  of  ordinary 
differential  equations.  In  order  to  apply  the  method  of  integral  relations,  the 
system  of  partial  differention  equations  is  integrated  across  the  layer,  after 
multiplication  by  appropriate  weighting  functions .  Hence  the  partial  derivatives 
with  respect  to  one  variable  are  eliminated,  and  in  the  resulting  ordinary  differ¬ 
ential  equations  there  appear  only  integrals  of  functions  of  die  dependent  variables . 
These  functions  are  approximated  by  suitable  interpolation  formulas . 

2.  Integral  Relations 

35 

To  obtain  the  integral  relations  for  die  system  of  transformed  equations  , 
equation  (V-4)  is  multiplied  by  a  weighting  function  F(U)  which  approaches  zero 
sufficiendy  rapidly  as  Y"  ®  ,  equation  (V-5)  is  multiplied  by  .  and  the 

resulting  equations  are  added: 


aK0F(0>]  +  sHwf<u>' 


%  dF(U)  2  dF(U) 
qb  dU  'U  dU  ay2 

(V-6) 


After  integrating  the  above  equation  with  respect  to  Y  from  0  to  ®,  one  obtains 


the  integral  relation 
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_d 
dX  J 


U  F(U)dY  = 


%  fdFQJ)  2  dF(U)  MJ_  I 
qb  JQ  dU  '  7  dU  Y  |  Y  =  0 


-J 


“  d2  F(U) 
0  dU2 


dY 


(V-7) 


This  equation  is  rewritten  in  terms  of  the  quantities 


■  (It)" 


and  0. 


0  \IY) 


-1 


Y  =  0 


(V-8) 


Changing  the  variable  of  integration  from  Y  to  U  then  gives 


_d 

dX 


r1  i  <i2f(u) 

J0  6  dU2 


(V-9) 


It  is  clear  that  0  approaches  infinity  as  U  -•  1.  From  the  Blasius  solution  for  a 

flat-plate  boundary  layer  with  no  pressure  gradient,  evaluated  as  the  similarity 
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variable  approaches  infinity  ,  one  obtains  a  singularity  of  the  form 


constant 


1 

<1-U) 


log  (1-U)  ] 


as  Y  -  ® 


In  most  previous  applications  of  the  method,  this  form  is  simplified  by  omitting 


the  logarithmic  factor  and  so,  for  U  -»  1,  0  is  approximated  by 


0  =  Ofd-U)'1  ]  as  U  ■*  1 

If  the  same  behavior  fir  assumed  in  the  present  case,  then  an  obvious  simple 
choice  for  the  weighting  functions  is 

fk  i<U)  =  (1-U)1  (V-10) 

As  in  other  applications,  the  functions  ©and  1/0 in  the  Kth  approximation  are 
represented  by  the  following  expressions 

e  =  <a0  +  AjU  +  . . .  +  ak1  uk*!)  (V-ll) 

0’1  =  (1-U)(B0+B1U+  ...  +  BK  1UK'1)  (V-12) 

The  coefficients  Aq,  A^,  . . B^,  B^  ...  are  functions  of  X  and  are  related  by 
the  condition  that  for 


where  J  =  0,  1,  2 . K  -  1,  the  representation  of  0  and  0  *  should  be  con¬ 


sistent: 
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A1  v 


\-iujK'1)Mb0^1uj+...+  Bk.i^-1)-1 


(V-13) 


When  K  =  1  the  approximate  expressions  for  0  and  0  are 


0 


e0(X) 

l-U 


01  =  ©qWU-u) 


(V-14) 


if 


F^U)  =  (l-U) 


dF:(U) 

dU 


d2F  (U) 

- =  0  (V-15) 

dU 


then  equation  (V-9)  becomes 


1 

0 


U0QdU 


qb  r1 

~Jo(i+u)0odu  + 


0. 


The  resulting  approximate  differential  equation  is 


dX 


3  — 
% 


6°  =  eo 


(V-16) 


which  may  be  integrated  directly.  The  solution  can  be  written  as 


(V-17) 
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Substituting  equations  (V-14)  and  (V-17)  into  equation  (V-8)  and  rearranging, 
one  obtains 


(V-18) 


Integrating  the  above  equation  with  respect  to  Y  for  fixed  X,  there  results 


(V-19) 


From  the  solution  of  the  inviscid  core  flow  the  velocity  component  in 
the  0  direction  is  and  Is  given  by  equation  (IV-30).  In  nondimensional  form 
it  can  be  written  as: 


p  =  2r0  (02  -0)  +  4 
q^_  dr  o  0 


g  |'exp[ -(2n+l)(l-r)7r/0o  ] 
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(2n+l)47r4J 


exp[  -  (2n+l)(l -r)7r/0o  ] 


sin(2n+l)  tt  0 


_  w 

4T 
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23.90 


(2n  tt)  J 


exp  [  -2n(l-r)*70  ] 


sin  2n  tt  0 


(V-20) 


The  velocity  of  die  inviscid  core  flow  along  die  base  can  be  evaluated  by  setting 
—2  2 

r=  1-0  (1-0  )/2,  and  is  given  by  the  expression 
0 
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qb  = 


20  (  62  -  0)  +  - - — -  sin(2n+l)7r' 

7 r  0l(2n+l)2. 


8  V 

+  2  L 
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7.40  360 
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7T  0 


i  /o_  ,v2  2  '  ..  , .4  4  2(2n+l) 

L  (2n+ 1)  7T  (2n+  1)  tt 


sin(2n+l)7r  G 
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*  l 


23.90  1 
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L(2n  jr) 
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( V-21) 


For  0  sufficiently  small,  X  can  be  approximated  by  r  0.  Since  also  r  0  , 

o  rr  3  o  o  o  r 


( V-22) 


Since  0  is  small,  Y  *  r  -  r  .  Then 
o  o 


Y  = 


/rT 

q  ^  j*  _ 


/Re 


qb(1-r) 


(V-  23) 


From  equations  (V-19)  through  (V-23),  the  solution  for  U  can  be  written  as: 


U  = 


1  -  exp 


(V-24) 
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When  K  =  2  the  approximate  expressions  for  ©  and  © 


-] 


©  = 


(1 


HT)  |>o<1-2U>+eiu’ 


If 


e'1  =  (l-U)[^-(l-2U)+|-  4u] 


dF  (U) 

F1  =  (1-U)'  ~ 


,  dF(U) 

f2  =  <1-U>  -  i£r 


d2F  (U) 

■i.  —4 —  =  o 


dU 


=  -2(1-U), 


d2F2(U) 


then  equation  (V-9)  becomes 


X  u[e0<i-2U)  +  eJo]du 


■  ‘  qT  J *  <1+U>  [  e0(1  "2U>  +  ®1 U  ]  dU  +  5" 


b  0 


u<i-u)[e0(i-2U)  +  etu]du 


2%  -1 


=  -—  J  (1-U2)  [e0(l-2U)  +  ©:  u]  dU 


b  0 


dU 


0  0 


0  1 


are 


(V-25a) 


( V-25b) 


(V-26a) 


(V-26b) 


(V-27) 


(V-28) 
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The  resulting  approximate  system  of  differential  equations  is 


d®° 

dX 


.2/96  x  70)3.32 

%  <9V79i>+e0  e, 


de, 

dX 


( V-29) 


(V-30) 


Since  and  $  can  be  obtained  from  the  solution  for  the  inviscid  core  flow 
evaluated  at  the  base  of  the  wedge,  the  above  systems  of  ordinary  differential 
equations  can  be  integrated  numerically  to  give  eQ(X)  and  G^X),  provided  that 
the  initial  conditions  can  be  obtained.  Since  the  first  of  equations  (V-25)  can 
be  rewritten  as 


dY*[<V2V(nc)+eo(iV)]dU  <v-31> 

then  for  X  fixed  the  above  equation  can  be  integrated  to  give 


Y  =  -  ( Ql  -  eQ)  fn(l  -U)  -  (  6:  -  29q)U  ( V-32) 

In  order  to  provide  initial  conditions  for  the  systems  of  approximate 
equations,  it  is  necessary  to  examine  the  behavior  of  these  equations  near  the 
wake  center  line.  For  0  sufficiently  small  the  velocity  of  the  inviscid  core  flow 
along  the  base,  given  by  equation  (V-21),  possesses  a  power -series  expansion  in 
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0  and  has  the  following  form: 

qb  -  c,e  +  c2e2+ . 

s  0  for  0  very  small 

Hence  from  equation  (V-22) 

X  a  iCj  02;  a  J2C^X  (V-33) 

• 

—  ®  for  8  very  small  (V-34) 

qb  2X 

For  this  value  of  q^/q^,  the  approximate  system  of  differential  equations,  which 
possesses  a  singularity  at  X  =  0,  can  be  shown  to  yield  an  exact  solution  of  the 
form 


=  Dj^/x"  for  X  «  l 


(V-35) 


The  constants  Dj  are  determined  by  the  solutions  of  the  algebraic  equations  which 
are  obtained  by  substituting  equation  (V-35)  into  the  original  approximate  system 
of  differential  equat  ons.  For  K  =  i,  then 
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e0  ■  Do7*" 


DQ  =  1  ;  ©0  =  Jx  for  X  «  1 


For  K  =  2,  then 


90  =  Doyx  :  e[  =  Dl/x 


10V  -  £0  -  ^ 


.n  4.  7n  -  —  M 
4D0  +  7D1  DQ  D1 


Solving  the  above  system  of  algebraic  equations,  one  obtains 


DQ  =  1.147  ;  Dj  =  1.822 


eQ  =  1 . 147  \[x  ;  ©1  =  1.822  Jx  for  X«1 


(V-36) 

(V-37) 

(V-38) 

(V-39) 


The  solution  given  by  equation  (V-39)  provides  the  initial  values  for  the  quantities 
©0  and  Gj,  if  evaluated  at  very  small  X.  Numerical  integration  of  the  systems 
of  ordinary  differential  equations  can  now  be  carried  out. 


VI 


LAMINAR  FREE  SHEAR  LAYER 
1  General  Description  and  Formulation 

The  simplified  model  of  the  shear -layer  flow,  which  retains  the  main 
physical  features  of  the  exact  flow,  is  described  here.  The  inviscid  supersonic 
flow  expands  at  the  trailing  edge  of  the  wedge  through  a  centered  Prandtl  -Meyer 
expansion  fan  to  the  base  pressure.  The  boundary-layer  flow  along  the  wedge 
surface  separates  at  the  trailing  edge  of  the  wedge  and  expands  isen tropically 
to  the  base  pressure  in  a  very  short  distance.  Subsequently  it  mixes  under 
nearly  constant  pressurewith  the  relatively  low-speed  recirculating  flow.  The 
expanded  boundary -layer  velocity  profile  forms  the  upper  part  of  the  initial  velo¬ 
city  profile  of  the  shear  layer.  The  remaining  part  of  the  initial  velocity  profile 
is  given  approximately  by  the  similarity  profile  obtained  by  assuming  a  mixing 
between  a  uniform  stream,  with  velocity  equal  to  the  velocity  at  the  inner  edge 
of  the  expanded  boundary  layer,  and  a  stagnant  region.  The  velocity  in  this 
portion  of  the  profile  is  evaluated  at  a  very  small  distance  from  the  trailing  edge 
of  the  wedge.  Near  the  dividing  streamline,  which  separates  the  recirculating 
flow  from  the  expanded  boundary  layer,  the  velocity  gradient  is  considerably 
larger  than  the  velocity  gradient  in  the  remaining  part  of  the  profile,  over  the 
entire  length  of  the  shear  layer.  The  thickness  of  the  shear  layer  is  of  higher 
order  than  r  ,  and  the  curvature  of  the  dividing  streamline  is  small  except  near 
the  stagnation  points . 
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From  the  description  of  the  simplified  shear -layer  model,  the  reie - 
tive  orders  of  magnitude  of  individual  terms  in  the  Navier -Stokes  equations  are 
estimated  in  the  same  manner  as  in  the  derivation  of  the  boundary -layer  equa¬ 
tions .  For  the  two-dimensional  case  the  resulting  differential  equations  are 
identical  to  the  flat -plate  boundary -layer  equations  with  pressure  gradient: 


Ml  .  Ml  =  o 

Sx  dy 


(VI-1) 


—  du  —  a  u 

pu  -=  +  pv  — =- 

ax  ay 


—  ah  —  ah 

p  u  —  +  p  v  r-^ 

ax  ay 


- 

dx 


(VI-2) 


(VI-3) 


where  x  is  measured  along  the  dividing  streamline  and  y  is  measured  normal  to 

the  dividing  streamline .  The  dividing  streamline  is  really  slightly  curved,  but 

the  curvature  affects  only  the  higher  approximations  to  the  equations  for  the 

mixing  region.  To  simplify  the  problem  further,  it  is  assumed  that  op  is  con- 

c  p 

stant  across  the  shear  layer  and  Pr  =  -j»-  is  equal  to  unity.  The  stream 
function  ip  is  defined  in  the  usual  manner  by 


p  Ml 

Me  ay 


and  a  modified  Howarth  transformation  is  introduced  as  follows: 
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i  - 
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x  ~  e  f  ~ 

=  r  :  *  =  -r  J„ p  dy 


o  0 


When  this  transformation  is  applied  to  equations  (VI-1),  (VI-2)  and  (VI-3), 
there  results  a  system  of  equations  similar  in  form  to  the  incompressible 
boundary -layer  equations: 


+  h  =  o 
a*  +  a?  0 


(VI-4) 


~  IV  IV 

~  flu  ,  ~  Ju  1  dp  flu 

u,sr  v»r 


(VI-5) 


(VI-6) 


Assuming  that  H  =  constant  in  the  initial  profile,  a  solution  of  equation  (VI-6) 
is 


~2 

~  ~  u 

h  =  H  -T 


(VI-7) 


which  reduces  equation  (VI-6)  to  (VI -5) .  This  result  implies  that  die  energy 


•  r 


I 
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equation  is  automatically  satisfied  once  the  momentum  equation  is  satified. 
The  density  in  the  shear  layer  can  be  written  as: 


Since  has  been  defined  as  the  pressure  immediately  after  the  corner  expan¬ 
sion  has  been  completed,  the  Bernoulli  equation  for  the  boundary  0  =  0^  of  the 
inviscid  core  flow  can  be  written  as 


Pe  =  P  +  i«r\ls2 


where  qg  is  the  nondimens ional  velocity  of  the  inviscid  core  flow  evaluated  at 
0  =  0q  .  Hence  the  density  in  the  shear  layer  can  be  rewritten  as 


P 


r\2> 

2  ~ 

(y-l)M‘(H  -  ~) 


(VI-8) 


The  momentum  equation  (VI -5)  is  therefore  independent  of  the  energy  equation 
(VI-6). 

33  34 

An  integral  method  of  the  Karman -Pohlhausen  ’  type  is  applied  to 


analyze  the  shear  layer,  because  it  is  simple  and  permits  the  treatment  of  the 
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case  when  the  shear  layer  is  not  thin  and  has  non -zero  inner -edge  velocity. 

This  method  involves  approximation  by  polynomials  with  undetermined  coeffi¬ 
cients  for  the  velocity  profiles.  The  boundary  conditions,  continuity  conditions, 
integrals  of  the  momentum  equation,  momentum  equation  along  the  dividing 
streamline,  and  momentum  equation  along  the  inner  edge  of  the  outer  shear 
layer  yield  a  sufficient  number  of  conditions  for  the  complete  determination  of 
the  profile  coefficients  and  thicknesses. 

2.  Velocity  Profile 

The  shear  layer  will  be  divided  into  an  outer  shear  layer  and  a  shear 
sublayer .  The  outer  shear  layer  is  the  portion  of  the  shear  layer  where  the 
kinetic  energy  of  the  fluid  is  large  and  the  velocity  gradient  is  not  too  large . 

The  initial  velocity  profile  is  determined  by  the  boundary -layer  velocity  pro¬ 
file  after  the  corner  expansion.  This  portion  of  the  shear  layer  grows  into  the 
inviscid  high-speed  outer  flow,  and  flows  downstream  of  the  rear  stagnation 
point  to  consitute  a  part  of  the  far  wake .  The  shear  sublayer  i;  the  portion 
of  the  shear  layer  with  larger  velocity  gradient,  in  the  vicinity  of  the  dividing 
streamline.  The  sublayer  flow  above  the  dividing  streamline  diffuses  into  the 
outer  shear  layer  and  flows  downstream  of  the  rear  stagnation  point  forming  a 
part  of  the  far  wake.  The  sublayer  flow  below  the  dividing  streamline  grows  in¬ 
to  the  relatively  low -speed  inviscid  recirculating  flow,  and  in  the  vicinity  of  the 
rear  stagnation  point  it  turns  back  toward  the  base  of  the  wedge .  The  initial 
velocity  profile  for  the  shear  sublayer  is  given  by  the  similarity  profile  resulting 
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from  the  mixing  of  a  uniform  flow  with  fluid  at  rest. 

Obviously  the  velocity  profile  has  different  character  in  the  different 
parts  of  the  shear  layer  and  the  dominant  factors  influencing  the  change  of  veloc¬ 
ity  profile  are  not  the  same.  Hence  it  is  difficult  to  describe  the  shear -layer 

velocity  profile  as  a  whole  by  a  single  expression .  In  the  study  of  mixing  be- 

36 

tween  a  separated  boundary  layer  flow  and  a  fluid  at  rest,  Kubota  and  Dewey 
represent  the  shear -layer  velocity  profile  by  two  polynomials,  one  used  above 
and  one  below  the  dividing  streamline,  because  the  initial  velocity  profile  is 
noticeably  different  from  the  subsequent  profiles.  Similarly  the  shear  layer 
discussed  here  will  be  divided  into  three  layers,  two  above  and  one  below  the 
dividing  streamline.  They  are  the  outer  shear  layer,  the  upper  shear  sublayer, 
and  the  lower  shear  sublayer  .  Each  layer  is  represented  by  a  polynomial  with 
undetermined  coefficients  as  functions  of  the  coordinate  measured  along  the 
dividing  streamline. 

Since  the  velocity  profile  for  the  shear -layer  flow  always  possesses  an 
inflection  point,  a  polynomial  of  at  least  third  degree  must  be  employed  to  repre¬ 
sent  the  portion  of  the  velocity  profile  where  the  inflection  point  is  most  likely  to 
37 

occur  .  It  is  known  that  there  is  no  point  of  inflection  in  the  initial  velocity  pro¬ 
file  for  the  outer  shear  layer.  Since  the  small  pressure  variation  induced  by  the 
recirculating  flow  and  the  relatively  small  viscous  force  have  little  effect  on  the 
outer  shear  layer  flow  where  the  kinetic  energy  is  very  large,  the  point  of  in¬ 
flection  will  occur  in  the  shear  sublayer . 
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For  simplicity  the  outer  shear  layer  may  be  approximately  represented 
by  a  second -degree  polynomial .  To  achieve  sufficient  accuracy  where  the  veloc¬ 
ity  gradients  are  large,  the  shear  sublayers  are  represented  by  fourth -degree 
polynomials.  The  velocity  profile  for  the  outer  shear  layer  is  expressed  as 

u  =  g(C)  =  aQ  +  aAC  +  a2?2  (VI-9) 

where  X  =  ^  and  ^  is  the  thickness  of  the  outer  shear  layer.  It  is  con- 

venient  to  introduce  the  quantity  Cq  because  of  the  small  displacement  effect  of 
the  shear  sublayer.  It  will  be  shown  in  section  5  that  ^  should  be  taken  to  equal 
6^/6  ,  where  6*  is  the  displacement  thickness  of  the  upper  shear  sublayer 
measured  in  terms  of  the  y  coordinate.  Hie  upper  shear  sublayer  is  repre¬ 
sented  by 


m  (V))  =  Lbn\+) 


(VI-10) 


where  =  -y-  and  6^  is  the  thickness  of  the  upper  shear  sublayer.  The 
lower  shear  sublayer  is  represented  by 


gtO)  =  ®(-)  (*(-))  =  ACn  V 

0 


(VI-11) 


where  ^  ^  and  6^  ^  is  the  thickness  of  the  lower  shear  sublayer. 
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The  profile  coefficients  aQ,  a 2>  b^,  ,  b^,  b^,  b  ,  cQ,  Cj,  c ^  c^, 

c^,  and  the  layer  thicknesses  5  ,  6^,  6^  ^  are  all  functions  of  x  .  They  are 
determined  bv  the  boundary  conditions  at  the  inner  and  outer  edges  of  the  shear 
layer,  the  continuity  conditions  between  the  layers,  and  appropriate  differential 
equations . 

Since  at  the  outer  and  inner  edges  of  the  shear  layer  the  velocity  must 
be  continuous  and  the  first  derivative  of  velocity  with  respect  to  y  must  be  zero 
(in  the  approximation  of  boundary -layer  theory),  the  boundary  conditions  are 


g(D  =  1 


(VI -12a) 


(VI-12b) 


and  the  boundary  conditions  at  y  =  '  ^  ^  are 


g(-)('°  ‘  (ffe  )(§{0) 


(VI -13a) 


*g(.) 

(-D  =  o 

V) 


( VI - 1 3b) 


In  order  to  have  a  smoother  fit  at  the  inner  edge  of  the  shear  layer  the  second 


derivative  of  the  velocity  with  respect  to  y  is  also  required  to  be  zero.  Hence 
an  additional  boundary  condition  at  y  =  ^  ^  is 
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a  g 


(-) 


~  9 

ac 

V) 


=  0 


(VI -13c) 


For  a  reasonably  smooth  fit  between  the  layers,  the  velocity  and  its  first  and 
second  derivatives  with  respect  to  y  must  be  continuous.  The  continuity  con- 

/s*# 

ditions  at  y  =0  are 


(VI-14a) 


(VI-14b) 


(VI-14c) 


It  does  not  seem  evident  whether  or  not  the  shear  sublayer  will  remain 

negligibly  thin  compared  with  the  outer  shear  layer  for  the  entire  length  of  the 

shear  layer.  Therefore  the  conditions  at  the  outer  edge  of  the  sublayer  are  not 

imposed  at  C  =  0  but  at  a  presumably  rather  small  positive  value  of  C  .  As 

will  be  shown  later,  through  the  requirement  of  equalization  in  pressure  across 

i 

x 

the  shear  layer  to  order  Re^  ,  the  dividing  streamline  must  be  displaced  rela¬ 
tive  to  the  outer  shear  layer  by  an  amount  equal  to  the  displacement  thickness 

of  the  upper  shear  sublayer.  Therefore  the  matching  point  between  the  upper 

V  ^(+)  V  ^(+) 

shear  sublayer  and  the  outer  shear  layer  is  not  at  C  =  -r  but  at  C  =  k  -*»— 7 , 
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where 


K 


As  noted  previously,  the  coordinate  y  is  measured  from  the  dividing  stream- 
line  but  C  is  defined  by  C  =  tp  -  t_ .  It  will  be  shown  in  section  5  that 

-j-O  0  0 

o  ~  ^ 

can  be  related  to  k  by?Q  =  (1  -  k)  =  . 


Since  the  inner  edge  of  the  outer  shear  layer  is  not  really  at  y  =0  and 
the  upper  shear  sublayer  is  assumed  rather  thin,  Taylor  expansions  about 
y  =  0  for  the  outer  shear  layer  may  be  used  to  provide  approximate  conditions 
at  the  matching*point.  Therefore  the  continuity  conditions  between  the  upper 
shear  sublayer  and  the  outer  shear  layer  can  be  approximated  by 


g(0)  +  k 


(VI- 15a) 


~2~ 

ML  <o)| 

a?2  J 

(VI-15b) 


(VI-15c) 


From  the  assumed  form  of  the  velocity  profile  given  in  this  section,  it  is 
clear  that  an  essential  feature  of  the  shear -layer  analysis  adopted  here  is  the 
division  of  the  shear  layer  into  three  parts .  Since  separate  quartic  profiles  are 
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used  to  represent  the  velocity  profiles  in  the  upper  and  lower  shear  sublayer, 
and  a  quadratic  profile  is  employed  to  express  the  outer  shear  layer,  there 
are  altogether  sixteen  unknown  profile  coefficients  and  thicknesses .  The  boun¬ 
dary  and  continuity  conditions  supply  eleven  relations  among  them.  However, 
since  the  total  number  of  relations  must  be  equal  to  the  number  of  unknown 
quantities  appearing  in  the  shear -layer  velocity  profile,  five  more  relations 
are  needed  so  that  the  velocity  profile  may  be  determined  uniquely. 

3 .  Momentum  Integral  Method 

An  integral  method  of  the  K^rm^n-Pohlhausen  type  is  employed  to  furnish 
some  of  the  remaining  relations  required  for  determining  the  unknown  coeffici¬ 
ents  in  the  assumed  representations  for  the  velocity  profile.  Three  relations 
may  be  obtained  by  taking  the  zeroth  moment  of  the  momentum  equation  for  each 
of  the  three  layers,  i.e.,  by  integrating  the  momentum  equation  over  the  thick¬ 
ness  of  each  of  the  three  layers . 

In  the  derivation  of  an  integral  condition  for  the  ojter  shear  layer  y  is 

***  "»  V  ~  ~  ~  y  d6 

temporarily  replaced  by  y  +  CqO  and  v  is  replaced  by  v  +  u  Cq  .  When 
expressed  in  these  transformed  variables,  equation  (VI-5)  has  the  same  form 
as  before,  and  is  to  be  integrated  from  0  to  6  : 


At  the  inner  edge  of  the  outer  shear  layer  the  new  v  component  of  velocity  is 
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very  small,  as  shown  in  section  5,  and  is  assumed  to  be  zero.  Then  integra¬ 
tion  by  parts,  by  using  equation  (VI -4),  gives  the  following  result: 


(VI-17) 


Substituting  equation  (VI-17)  into  equation  (VI-16),  there  results 


P  ~2  ~  ** 

i  <u  'uHy 
0  x 


(VI-18) 


From  equations  (VI -9),  (VI-12b)  and  (VI-18)  the  zeroth  moment  of  the  momentum 
equation  for  the  outer  shear  layer  can  be  written  as 


By  using  equation  (VI-8),  p  =  p  (x),  and  the  Bernoulli's  equation  for  the  boundary 
0  =  0q  of  the  inviscid  core  flow,  equation  (VI -19)  is  reduced  to 


dq,  ~  ] 

I‘[B -4s2]ar 


1  -  ^  M2  q2 
2  r  ws 


(VI- 20) 


0 


For  the  upper  shear  sublayer,  equation  (VI-5)  is  integrated  from  0  to 


(+) 
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~  ~  ~  p  (+\  ~  ~ 

u  u~)dy  +  J  (vu~)dy 


y  y 


)  dy 


(VI-21) 


where  y  and  v  are  again  measured  relative  to  the  dividing  streamline.  Since 
v  =  0  along  the  dividing  streamline,  integration  by  parts,  using  equation  (VI -4), 
gives  the  result 


(VI-22) 

~ 

where  u^  is  the  velocity  at  the  outer  edge  of  the  upper  shear  sublayer  .  Substi 
tuting  equations  (VI -8)  and  (VI-22)  into  equation  (VI-21)  and  using  Bernoulli's 
equation  as  before,  then 


r’(+)  r,~2,  ~  ~  "I  ~  3u  I  (+) 

(u  r*  -u  u~  dy  =  -T7  1 
Jn  L  x  u  x  J  3y 


(y-l)M^  qo  dq: 


r  "s 


dST 


1  -  q2 
2  r  ^s 


J  (H  -  iu2)  dy  (VI-23) 


which  may  be  rewritten  as 


6 

(+)  , 

IS-  )2 

0 

\f°>/ . 

dy  + 


1  dg 
g(0)  d* 
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Uu  V  U  \  ~  _1_  /  _u_  \ 

mo)/\  g<o)/  J ay "  m  ay  \  m)  0 


ql  y 

_ dx 

1  -  ^  M2  q2  Jo 

2  r 


r<*>  r_H_ 

Jo  i  g  2<0)  Vv  J 


(VI -24) 


As  implied  by  the  definition  (VI -9),  g(0)  represents  the  approximate  velocity 
profile  for  the  outer  shear  ltyer  evaluated  at  Z  =  ^  =0,  and  is  a  function  of 
x  .  From  equations  (VI-9),  (VI-10),  and  (VI-24)  and  through  tie  application  of 
the  Leibnitz  rule,  the  zeroth  moment  of  the  momentum  equation  for  the  upper 
shear  sublayer  can  be  written  as 


d5r(.V)  J0®<+)  dV)j  '  W1}  d5c(  V)V(+)’ 

(*f0)  '^+)(1)^+)]d^( 


’(+)  V). 


<+>)  + 


(Y-DM2  q„  TXT 


1  ~ 


r  Ms  d5T  (+)  T  T  H  ,~2  l,y  , _ _ 

-O -  U^T  -**<+>  J‘,r<+>  (v,'25) 


1  *  9 

2  r 


0  g  (0) 


For  the  lower  shear  sublayer,  equation  (VI-5)  is  integrated  from 


0  B0 

***  ***  ***  r  ***> 

~  (u  u~)  dy  +  (v  u~)  dy 

J-6  x  y 

V)  (-) 

=  -  idy  +  ( IT - )  d y 

d5f  J  y  V  3  Jr  y  y 

'  (-)  <-) 


(VI -26) 


to  0: 
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Again  since  v  =  0  along  the  dividing  streamline,  the  integration  by  parts,  by 
using  equation  (VI-4),  gives  the  following  result: 


Substituting  equations  (VI-8)  and  (VI-27)  into  equation  (VI-26)  and  using  Bernoulli's 
equation  as  before, 


•o 

V 

(“?  x  ‘(r)  %(“)? 

t 

i 

e 

r 

a_u 

ay 


0 

-y  + 

V) 


2  dqs 
(y  ~!)  Mr  Qg 

1  -  ^M2q2 
2  r^s 


iu  ^  dy 

(VI-28) 


which  may  be  rewritten  as 


From  equations  (VI-9),  (VT-11),  (VI-13a,b)  and  (VI-29)  and  through  the  applica- 
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tion  of  the  Leibnitz  rule,  the  zeroth  moment  of  the  momentum  equation  for  the 
lower  shear  sublayer  can  be  written  as 


dx(V)I.,  g(-)dV)  )  "(  Ug \%0))  d*( V)  (j  g(-)dV)  ) 

=  -  ‘^o) 


■(-) 


2  d%  ~ 

(,-DMr  U»,., 

1  -  *  M2  q2 
2  r  Hs 


0/  ~ 
f— 

l  ~2  2 

J-Ag  (0) 


+  (  1  »*-> 


0 

-1 


(VI -30) 


Each  of  the  preceding  integral  conditions  requires  an  overall  momentum 
balance  in  one  of  the  three  parts  of  the  shear  layer .  Two  more  conditions  are 
obtained  by  requiring  that  the  momentum  equation  also  be  satisfied  at  particular 
locations  in  the  shear  layer.  Along  the  dividing  streamline  the  momentum  equa¬ 
tion  (VI-5),  with  u  (x,  y)  replaced  by  g(0)g^(0)  according  to  equation  (VI-10) 


becomes 


dg,^(0> 


_L  d  g(0)  \ 

^  rv\  j(V  I  ' 

2  *s 


inP  -W°>(p)  W) 


dq 

'-XA% 


(^+)(0)e<0,6(2+)  ) 


/  H  J_  ~2  \  _1 

(;2<o)  2g<+>(  )  v/0) 


3  2  a 


(+) 


<> 


(VI-31) 


(0) 


The  flow  along  the  inner  edge  C  =  0  of  the  outer  shear  layer  is  also  required 
to  satisfy  the  momentum  equation.  Then  from  equation  (VI-5)  and  (VI -9)  one 
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obtains 


d_im> ,  j  i!i(0)  + 
d*  ?(o>y2  «2 


2  dqs 

W 

1  -  *M2q2 
2  r 


[  H  •  is2(°)J  §5) 


(VI -32) 

In  addition  to  the  eleven  relations  already  obtained  from  the  boundary 
and  continuity  conditions,  die  ordinary  differential  equations  (VI-20),  (VI-25), 

(VI-30),  (VI-31)  and  (VI-32)  are  used  to  supply  the  required  independent  rela- 

# 

tions  for  the  complete  determination  of  the  velocity  profile  in  the  shear  layer. 
However,  it  is  also  necessary  to  provide  the  initial  velocity  profile  in  the  shear 
layer. 


4  .  Trailing-edge  Expansion  and  Initial  Development  of  Shear  Layer 


The  expansion  process  at  die  trailing  edge  of  the  wedge  is  assumed  isen- 
tropic  and  the  total  enthalpy  is  regarded  as  constant  throughout  ‘he  boundary 
layer.  It  can  be  shown  that  the  velocity  after  corner  expansion  is 


(VI -33) 


where  f  is  the  Blasius  solution  for  an  incompressible  flat -plate  boundary  layer 
and  the  subscript  c  refers  to  the  condition  at  the  trailing  edge  of  the  wedge  and 
at  the  outer  edge  of  the  boundary  layer. 

If  the  length  of  the  wedge  surface  is  denoted  by  Lc>  then  the  Blasius 
similarity  variable  evaluated  a  very  short  distance  upstream  from  the  trailing 
edge  of  the  wedge  may  be  represented  by 
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(VI -34) 


where  y  is  the  y  coordinate  before  the  expansion,  just  upstream  from  the  trail¬ 
ing  edge  of  the  wedge.  Since  there  is  no  pressure  gradient  across  the  boundary 

~  2 


layer,  and  H 


the  relation  between  y  and  n  is 

c  c 


H  -  i 


(VI -35) 


where  f'  =  f qC T) c )  -  Since  the  total  mass  in  the  boundary  layer  immediately 
before  and  after  the  comer  expansion  must  remain  the  same,  the  relation  be¬ 
tween  y  and  t)  is 
e  'c 


(VI -36) 


where  yg  is  the  y  coordinate  measured  from  the  dividing  streamline  after  the 
expansion  a  very  short  distance  downstream  from  the  trailing  edge  of  the  wedge. 

t 

Hence  the  relation  between  the  transformed  coordinate  and  the  similarity  vari¬ 
able  is 


/  2Re  /  L  p  u  \  f' 
~  /  _ e  \  (  c  c  c  c 

y  \Re  /  \  r  pu/TT 

'  /-»  '  p  p 


dT?, 


( Vl-37) 


o  e  e 


where  Rec  =  u^  L  /Vc  .  From  equations  (Vl-33)  and  (Vl-37)  the  velocity  profile 
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after  the  trailing-edge  expansion  may  be  obtained  in  transformed  coordinates. 

The  boundary -layer  flow  immediately  after  the  trailing-edge  expansion 
mixes  with  the  flow  below  the  dividing  streamline  which  is  virtually  stagnant 
when  x  is  very  small,  i.e.,  very  close  to  the  wedge  base.  Hence  the  govern¬ 
ing  differential  equation  is  the  same  as  equation  (VI-5),  except  that  no  pressure  - 
gradient  term  is  present,  and  can  be  written  in  terms  of  the  stream  function  as 

2~  2~  3~ 

|l  M  =  ijL  (VI-38) 

Sv  ox3y  Sx  , ~ 2  ~3 

ay  ay 

For  the  outer  shear  layer  the  stream  function  is 


*l>  3  i^(y) 


where  ^(y)  =  '  ^en  aQ^’  a^0^’  311(1  ^(0)  are  °^tained 

by  fitting  a  polynomial  to  the  velocity  profile  after  the  corner  expansion.  For 
the  flow  near  the  dividing  streamline  the  viscous  stresses  are  very  large  and 
the  flow  is  described  approximately  as  the  mixing  of  a  uniform  flow  with  fluid 
at  rest.  Hence  it  is  assumed  that  the  solution  has  the  approximate  form 


$  ~  J 2aQ(0)x  fQ(T?) 


(VI-39) 
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n 


Fp 

Jit  y 


(VI-40) 


The  solution  is  obtained  from  the  following  differential  equation  and  boundary 
conditions 


f0<®)  1  :  0  frf°>  =  0 

Therefore  for  a  fixed  small  positive  x  and  for  t)  2  0, 

^♦)(r<+>)  “ 

Similarly  for  ^  s  0, 


(Vl-41) 

(VI-42) 


(VI-43) 


S->(r(->)  “  fo<?>  <v'-44) 

Initial  values  are  provided  for  the  profile  parameters  of  the  upper  and  lower 
shear  sublayer  by  fitting  polynomials  to  the  upper  and  lower  parts  of  f^(rj  )• 


S .  The  Location  of  Dividing  Streamline 


Since  the  velocity  in  the  inviscid  core  flow  is  not  exactly  zero,  the  pres- 

sure  is  not  exactly  constant,  and  therefore  the  shear  layer  is  not  exactly  straight. 

In  a  first  approximation  the  curvature  of  the  shear  layer  will  be  such  that  the 

pressure  in  the  external  flow  balances  the  pressure  in  the  recirculating  flow. 

The  velocity  in  the  inviscid  core  flow  is  •=-  =  0(M  )  and  so  the  nondimen sional 

2  e 

shear-layer  curvature  is  0(M^).  If  the  pressure  gradient  is  known  from  the 
solution  of  (IV),  the  shear -layer  curvature  can  be  calculated  to  this  order. 
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The  pressure  variation  across  the  shear  layer  must  be  zero  not  caly  to 

2  38 

order  but  also  to  higher  order  in  Reynolds  number.  Lock  ,  in  his  study 

of  the  mixing  of  parallel  streams,  notes  that  the  dividing  streamline  is  not  nec  - 

39 

essarily  straight.  Ting  shows  that  for  the  mixing  of  a  uniform  supersonic 
stream  with  a  parallel  uniform  subsonic  stream  the  dividing  streamline  is  dis¬ 
placed  toward  the  subsonic  portion  of  the  flow  by  an  amount  equal  to  the  dis  - 
placement  thickness  of  the  viscous  layer  above  the  dividing  streamline,  so  that 
the  pressure  is  kept  the  same  across  the  mixing  layer  up  to  an  appropriate  or¬ 
der  of  magnitude.  The  location  of  the  dividing  streamline  is  important  in  the 
study  of  the  shear  layer  in  the  wake,  because  it  has  an  influence  on  the  shear- 
stress  distribution  and  velocity  distribution. 

Since  the  velocity  distribution  in  the  outer  shear  layer  and  in  the  recir¬ 
culating  flow  is  very  complex,  it  seems  necessary  to  retain  the  approximation 
that  is  small.  It  is  assumed  that  the  outer  shear  layer  may  be  approximated 
by  a  uniform  shear  flow  at  supersonic  speed 

u  =  u  +  fiy 

where  is  the  vorticity  and  is  a  constant.  This  approximation  is  consistent 
with  previous  assumptions,  since  only  die  portion  of  the  outer  shear  layer  near 
y  =  0  needs  to  be  correctly  represented  here.  The  recirculating  flow  is  assum¬ 
ed  to  have  zero  velocity.  Since  the  recirculating  flow  will  actually  have  velo- 
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city  -  0(M  ),  then  p  =  p  [  1  +  0(M^)  3  in  the  recirculation  region,  and 


the  velocity  at  the  inner  edge  of  the  outer  shear  layer  will  actually  be 
_  _  2 

u  =  u  [  1  +  0(M  )  ]  .  The  nondimens ional  variables  for  the  uniform  shear 
u  r 

flow  are  defined  as 
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Or 


The  nondimensional  variables  for  the  recirculating  flow  are  defined  as 


x  y 
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r  7  r 

o  o 


l()  u  v 
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u  u 
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The  nondimensional  variables  for  the  mixing  layer  are  defined  as 
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+  x  +  V _ u  n  u 
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Following  the  same  procedure  as  for  the  mixing  of  two  parallel  streams 
39 

given  by  Ting  ,  the  flow  variables  may  be  expanded  as 


a  * 

V  'V- 

°0  +  ‘°1  +  • 

<v  <v 
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p  -pe 

~  "o  +  ‘»i  + 
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}  (VI-45) 


u  ~  Uq  +  € +  .  . 

2 
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\  (VI-46) 


>  (VI-47) 


where  C  =  Re  ^  ,  p  is  the  density  of  the  uniform  shear  flow,  and  p_ 
u  u  0 

=  0(M^). 


P0(x) 


i 


% 
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Substituting  equation  (VI -47)  into  the  Navier-Stokes  equations  and  collect¬ 
ing  coefficients  of  like  powers  of  <  gives  the  following  form  for  the  yt  com¬ 
ponent  of  the  momentum  equation: 

‘Pfyt  +  *  °<(2> 

The  matching  conditions  for  the  pressure  at  the  outer  and  the  inner  edge  of  the 
mixing  layer  give 

pt(x,  »)  =  {^(2,0) 

P|  (x  ,  -“)  =  Pj(x, 0) 

There  results  a  compatibility  condition 

^(3,0)  =  Pj(x,  0)  (VI-48) 

Equation  (VI-48)  states  that  the  pressure  difference  across  the  mixing  layer 
should  be  0(  €  ).  This  condition  is  used  at  a  later  stage  for  the  determining  of 
the  location  of  the  dividing  streamline  relative  to  Jie  outer  shear  layer. 

The  nondimens ional  Navier-Stokes  equations  for  incompressible  flow  may 
be  written  as 


■a  .  a  \  2 ,  222, 

T  Jx  x  ay  / 


(Vl-49) 
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From  equations  (VI-45)  and  (VI-49)  the  following  set  of  equations  can  be 
obtained: 


(*o*  ^  ■  $ok  =  0  <yl-50o) 

n  ' *o* ^ )  A +  ( *1}  A  •  *i3  * 0 

(VI -50b) 

Equation  (50a)  can  be  reduced  to 


2A  A 

V  Vq  =  N  =  constant 


(VI-51a) 


Since  the  following  boundary  conditions  must  be  satisfied 


$q(*'9)  =  $+  iNy2  upstream 


5>0<a,o)  =  0 


(VI-51b) 
( VI-51c) 


the  solution  is 


i  a.  a  N  a2 

*0(x-y>  =  y+  2y 


(VI-52) 


2 

As  noted  previously,  terms  CKM*  )  in  the  coefficient  of  9  have  been  neglected. 


Substituting  equation  (VI-46)  into  equation  (VI -49),  there  results  a  set  of 
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differential  equations 


(*0y  ty  '  *0*  #)A»-  0 


(VI -53a) 


r,  J. 


\St  '  ^Ox  ^  V  *1  +  (  *ly  to  •  *lx  ^  )V  ^0  =  0 

( VI- 53b) 

Since  ■=-  =  0(M  )  in  the  recirculation  region,  equation  (VI-5Sj)  can  be  written 
e 

as 


V  l i)Q  =  0(Mr) 


( Vl-54a) 


Since 


t^0(x,O)  =  0 

in  the  present  approximation  the  solution  will  be  estimated  as: 


(VI -54b) 


»J>0(x,  y)  =  CKMr)  (VI-55) 

Substituting  equation  (VI -4 7)  into  equation  (VI -49),  there  results  a  set  of 
differential  equations 


gyF (  *t ytytyt  +  Vxt  *  lytyt  ■  *lV  Xtyt  )  =  0 

(VI -56a) 

3yt  (  ^  2yt  yt  yt  +  ^ lxt  ^2yt  yt  *^lyt  *^2xtyt 
+  ^  2xt  ^lytyt  ‘  ^+2yt  ^Ixtyt  )  =  ° 


(VI  56b) 
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The  matching  condition  as  y  t  -*  ®  and  y  0  is 


0liyt(xt  «)  =  0Oy(x,  0) 


The  matching  condition  as  yt  -*  -  «  and  y-  0  is 


^+lyt  (x+ '  •  =  %y(X<0)  =  0 


Hence  the  following  asymptotic  relations  are  obtained 


$tj(xt,yt)  ~  ytjwx.0)  +  o(yt) 


as  yt  -  oo 

(VI -57a) 

as  yt  -»  -  oo 

(VI*57b) 

By  using  equation  (VI-57),  equation  (VI-56a)  can  be  reduced  to 


j/it^xt  ,yt)  ~  o(yt) 


(VI-58) 

Since  ^(x,  y)  is  a  function  of  y  alone,  equation  (VI-58)  can  be  reduced  to 


^lytytyt  +^lxt  ^lytyt  ‘  *  flyt  ^lxtyt  =  ° 


(VI-59) 
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Two  of  the  boundary  conditions  are 

<*t,  «)  =  1  (VI- 60a) 

0 tiy+  (set,  -«)  =  0  (VI -60b) 

Equations  (VI-59)  and  (VI-60)  have  a  solution  of  the  form 

=  rjt  = 

where  f(Tjt )  satisfies  the  differential  equation  and  boundary  conditions 


f +  f  f"  =0 

(VI-61a) 

f '(•>  =  1 

(VI -61b) 

f'(-«)  =  0 

(VI- 61c) 

Equation  \ VI-58)  actually  applies  to  the  mixing  layer  with  the  stream¬ 
line  =  0  slightly  displaced  relative  to  the  outer  shear  layer.  Let  xt  and  yt 
be  measured  along  and  normal  to  fp  ^  -  0  respectively.  Then,  since  the  curva¬ 
ture  of  the  line  «/)t  =  0  affects  only  the  higher  approximations  to  the  equation 
for  the  mixing  region,  equation  (VI-59)  remains  unchanged.  A  third  boundary 
condition  is  obtained  from  the  requirement  that  the  pressure  difference  across 
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the  mixLig  layer  be  of  order  *  .  Therefore,  the  velocity  profiles  in  the  first 
approximation  of  the  mixing  region  should  be  so  oriented  that  the  pressures  in¬ 
duced  in  the  inviscid  stream  by  the  effect  of  displacement  thickness  are  balanced 
across  the  mixing  layer.  Since  equations  (VI-61a,b,c)  are  invariant  under  the 
transformation 

T|t  -  Tjt  +  p t 

the  third  boundary  condition  may  be  assumed  as 


f<0)  =  0 


(VI-61d) 


The  solution  obtained  by  using  boundary  condition  f(0)  =  0  differs  from  those 
using  f(/3t)  =  0  only  in  the  location  of  the  line  ^  =  0.  The  solution  of  equa  ■ 

tions  (VI -61a, b,  c,d)  for  T)t  •  is 


f(Tjt)  ~  (|jt  +  +  exponential 


(VI-62) 


Where  /jt^  is  proportional  to  the  displacement  thickness  of  the  upper  part  of  the 
shear  sublayer  which  is  measured  from  =0.  Similarly  the  solution  for 
tj  t  -»  -  ®  is 


f(rjt)  ~  (rjt  +  /jtj)  +  exponential 


(VI-63) 
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where  /3'f.  is  proportional  to  the  displacement  thickness  of  the  lower  part  of  the 
shear  sublayer,  which  is  measured  from  =0. 

By  using  equations  (VI -53b),  (VI -54a)  and  (VI -63)  the  flow  below  the  mix¬ 
ing  layer  due  to  displacement  thickness  is  described  by 

V2*^  =  0  (VI -64a) 

j/j^x.O)  =  «/ 2x  for  x  >  0  (VI-64b) 

From  Bernoulli's  equation  p^  may  be  expressed  as 

Pj  =  -uQu1  (VI -65a) 

Since  € ,  <u,  and  €  v,  are  of  order  c  and  u„  is  of  order  M  ,  c  p.  will  be  esti- 
^111  0  r  *1 

mated  as 

cp:  =  0(Mr€) 

=  o(c)  (VI -65b) 

Therefore  one  gets 

p^x.O)  =  p^x.O)  =  ptj  (xt,yt)  =  0  (VI-66) 

if  terms  of  order  are  neglected  compared  to  unity. 

For  the  uniform  shear  flow,  the  next  approximation  is  governed  by  the 
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0 


(VI-70) 
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Hence  equation  (VI -68b)  reduces  to 

N0  =  0  (VI- 71) 

A 

3'nce  N  =  constant,  it  is  required  for  the  uniform  shear  flow  at  supersonic 
sDeed  that 


^($,0)  =  0^ (x, 0)  =  vtj (xt ,  co)  =  0 


(VI -72) 


This  relation  means  that  the  line  tjj  t  =  0  is  displaced  by  an  amount  equal  to  the 
displacement  thickness  of  the  upper  shear  sublayer.  In  other  words,  it  is  re¬ 
quired  that  there  be  no  effect  of  displacement  thickness  on  the  uniform  shear 
flow.  If  this  result  is  applied  to  the  problem  discussed  in  the  previous  sections, 
one  can  make  the  approximation  that  for  y  measured  from  the  dividing  stream¬ 
line,  the  inner  edge  of  the  outer  shear  layer  is  at  y  =  ?q6>  where  Tq  =  , 

and  the  v  component  of  velocity  is  negligibly  small  at  the  inner  edge  of  the  outer 
shear  layer . 

It  is  evident  that  the  approximation  (Vl-72)  is  not  highly  accurate,  be¬ 
cause  the  error  is  CKM^).  However,  the  only  purpose  here  is  to  obtain  an  esti¬ 
mate  for  k  ,  which  in  turn  will  be  used  to  estimate  the  effect  of  nonzero  sublayer 
thickness  on  the  outer  shear  layer.  Since  this  effect  is  expected  to  be  small,  it 


is  not  necessary  that  the  value  of  k  be  extremely  accurate. 


VII 


NUMERICAL  SOLUTION 
1 .  The  Derivation  of  the  Equations  for  Computation 

In  the  discussion  of  the  conservation  of  angular  momentum  for  the  en¬ 
tire  recirculating  flow,  it  was  observed  that  for  a  steady  motion  the  accelerating 
torque  from  the  shear  layer  should  always  be  equal  to  the  retarding  torque  from 
the  base  boundary  layer.  In  order  to  determine  these  torques  it  is  necessary  to 
find  the  distribution  of  the  shear  stress  both  along  the  base  of  the  wedge  and  a- 
long  the  dividing  streamline. 

Since  the  distribudon  of  the  shear  stress  along  the  base  of  the  wedge 
can  be  determined  by  the  solution  of  equations  (V-l)  and  (V-2)  for  the  base 
boundary  layer,  the  solution  of  the  system  of  approximating  ordinary  differ¬ 
ential  equations  (V-29)  and  (V-30)  is  needed.  For  fairly  good  accuracy  the 
velocity  of  the  inviscid  core  flow  evaluated  along  the  base  of  die  wedge  can  be 
approximated  by 
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because  the  series  in  q^  converge  quite  rapidly.  The  transformed  coordinate 

along  the  X  direction,  given  by  equation  (V-22),  can  be  approximated  by 

0 

X  =  a  d0.  Hence 
0  D 


%  69 


(VII -2) 


Equations  (V-29)  and  (V-30)  can  be  rewritten  as 


(VII -3) 


(VII-4) 


These  equations  are  ordinary  differential  equations  with  0  as  the  independent  vari- 

2 

able.  The  initial  conditions  for  ©q  and  ©^  can  be  obtained  by  substituting  X^  0 
into  equation  (V -39).  Then,  for  0  fixed  but  small,  the  following  approximation  is 
obtained 


©Q  a  0.815  ( VII -5) 

e  a  1 . 292  «/C^  0  (VII-6) 

where  may  be  obtained  from  the  Taylor  expansion  of  q^  about  0=0,  and 
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Hence  the  Runge-Kutta  method  can  be  used  to  provide  solutions  for  ©Q  and  0^ . 


Since 


/Re" 


u  =  Ji 


and  0 


au 

aY  Y  =  0 


the  shear  stress  along  the  base  of  the  wedge,  T  ,  can  be  written  as 


*  (VII-7) 


the  local  skin -friction  coefficient  along  the  base  of  die  wedge,  C_„  has  the  form 
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(vn-8) 


The  distribution  of  the  shear  stress  along  the  dividing  streamline  can  be 
obtained  from  the  solution  of  the  shear-layer  equations  (VI-1),  (VI-2)  and  (VT-3). 
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An  integral  method  of  the  Karman-Pohlhausen  type  is  used  to  reduce  this  set  of 
partial  differential  equations  to  a  set  of  approximating  ordinary  differential  equa¬ 
tions  .  The  velocity  profiles  are  approximated  by  polynomials  with  undetermined 
coefficients,  given  by  equations  (VI -9),  (VI-10),  and  (VI-11).  Since 
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the  distribution  of  the  shear -stress  along  the  dividing  streamline,  T  can  be 
written  as 
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The  local  skin  friction  coefficient  along  the  dividing  streamline, c^,  has  the  form: 
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(VII -10) 


where  a^,  b^,  and  6^  are  given  by  the  solutions  of  the  system  of  equations 
(VI-20),  (VI-25),  (71-30),  (VI-31)  and  (VI-32)  with  the  boundary  and  continuity 
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conditions  (VI -12)  through  (VI- 15). 

Since,  in  the  continuity  condition  between  the  outer  shear  layer  and  the 
upper  shear  layer,  the  location  of  the  outer  edge  of  the  shear  sublayer  is  need¬ 
ed,  the  location  of  the  dividing  streamline  must  be  determined.  It  has  been  shown 
that  the  dividing  streamline  is  displaced  by  an  amount  equal  to  the  displacement 
thickness  of  the  upper  shear  sublayer.  This  satisfies  the  requirement  of  zero 
pressure  gradient  across  the  shear  layer  to  order  (Re^). 

The  displacement  thickness  for  the  upper  shear  sublayer,  is  ex¬ 


pressed  as 
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where  u^  =  u^g(0)  is  the  velocity  at  the  outer  edge  of  the  upper  shear, u^  =  u  g(tcA) 

V  l~ 

is  the  velocity  at  the  outer  edge  of  the  upper  shear  sublayer,  A  =  (+)  /  6 

The  outer  edge  of  the  upper  shear  sublayer  is  at 


C  =  k  A 


( VII  - 1 2) 
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/0\{T*  /  6* 

where  k  =  (  l  -  \  Since  can  be  shown  to  be  nearly  constant,  -pif 

\  6<+)/  *o>  «r+) 

may  be  assumed  negligibly  small  in  the  computation.  This  will  simplify  the  compu¬ 
tation  and  is  within  the  accuracy  needed.  Hence  equations  (VI-15a,b,c)  can  be  re¬ 


written  in  the  following  form: 
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For  simplicity,  define 
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K  2  =  'K0A 


(VII- 15) 
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(VII -18) 


(VII- 19) 


1+A 


These  quantities  are  functions  of  an,  6,  6,  .,  6.  .  and  q  .  The  velocity  of  the 
inviscid  core  flow  evaluated  along  the  dividing  streamline  can  be  approximated 
by 
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Therefore,  based  on  the  boundary  and  continuity  conditions,  the  profile 
coefficient  can  be  expressed  in  terms  of  a^,  b^,  6,  6^  and  0^  : 
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where  6  =  ■ 

V) 

Substituting  equations  (VII -9),  (VII  *21)  and  (VII -22)  into  equations  (VI -20) 
and  (VI -31),  there  results 
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dq. 


(Y-l)Mr2qg 

where  k  =  - - — - -  .  Equations  (VI-32)  and  (VI-33)  are  two  differential 

1  -*M  q 
2  r  Ha 

equations  in  two  dependent  variables .  The  initial  value  for  6  and  a^  can  be  approxi 
mated  by  appropriate  curve  fitting  of  the  expanded  velocity  profile,  which  is  found 
by  assuming  an  isentropic  expansion  to  the  base  pressure  of  the  Blasius  profile 
evaluated  at  die  trailing  edge  of  the  wedge.  Ibis  expanded  velocity  profile  can  be 
obtained  from  die  solution  of  the  following  equations: 
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(VII -35) 


(VH-36) 


Equations  (VII-34),  (VII-35)  and  (VH-36)  are  obtained  from  the  Blasius  equation. 
Equations  (VII-37)  and  (VII -38)  are  obtained  from  equations  (VI-33)ihroigi(VI-37) 
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and  y  =— -  .  The  initial  conditions  for  the  system  of  differential  equations  (VII-34), 
c  lc 

(VII -35),  (VII -36)  and  (VII -3$)  are 

MO)  =  f(1)(0)  =  0 
c  c 

f  =  0.4696 
c 

y(0)  =  o 

Hence  the  Runge-Kutta  method  can  be  us  2d  to  provide  solutions  for  u  (rjc)  and 
y (tjc)-  A  plot  of  TT(y)  gives  the  expanded  velocity  profile  in  the  coordinate 
system  of  the  modified  Howarth  transformation.  For  x  =  0,  it  is  assumed 
that 

u<7>  a  1  -  (1  -  a0(0))(  1  -  jLf  (Vn-39) 

Therefore  a^O)  and  MO)  can  be  determined  approximately  from  the  plot  of  the 
expanded  velocity  profile.  The  system  of  ordinary  differential  equations  (VII-32) 
and  (VII-33)  can  be  solved  numerically  by  using  the  Runge-Kutta  metho^.  . 

Substituting  equations  (VI -9),  (VI-10)  and  (VI-11)  into  equations  (VI-25), 
(VI-30)  and  (VI-32),  there  results 


»  { V)  [  b0 +  Vl  +  5  Kb2 +  bl2  )  +  1  (  b0b3  +  bl  b2  )  + 


H2b0b4  +  b2  +  2blb3)+Kblb4  +  b2b3) 


4  (b32  +  2b2b4  )  +  ib3b4  +  9  b42l  } 


(I'ao>A  \  d 


-  115  - 


■  U^'m) 


2bO  +  2bObl  +  !(2bOb2  +  bl2) 


+  (b0b3  +  blb2  )  +  5  (2bOb4  +  b2  +  2blb3  ) 

-  I  (  blb4  +  b2b3  )  +  f  (  b32+  2b2b>  *b3b4  +  I b4 

^0|(V)[^-(VX^)!-,} 


4  3  4 


_d 

d2f 


'<+> 


,2 

:o 


■H-i 


itfXiihA'.-b*  4  4)1 


fl  i 

2  2 

f  b0  +  bObl  +  1 

(2b0b2+bl2) 

J 

lao 

L 

3  +  bl  b2  )  +H2bOb4  + 

b2b3)+?(b32+2b2b4) 

•*“*'11 

> 

(vn-40) 

:0C1  +  3  ( 

f  2  \ 

,  2c0c2 + ci ; 

+  C1C2) 

^  j  |  2  ^0^4  ^ 

C22+2C1C3  ) 

+  C2C3  ) 

+  i(c2+2c2 

C4  )  ?C3C4  +  iC42]} 

-  116  - 


--(r0  dF  «(+))  [2c0-2Vi  +  I(2C0C2+C12)-(C0C3+C1C2) 


+  5  (2C0C4  +  c2  +  2cic3)  '  3  (cic4  +  c2c3  )  ?(c3  +  2c2C4) 

•5C3C4+fC42]  '[(s^K)]  (C0'T  +  /‘-r+  f) 


+  +  '5  [  C0  '  C0C1  +  J  (  2t0C2+  c12  ) 

•I  (  C0C3  +  clc2  )  +  H2C0C4  +  c2  +  C2  +  2C1  c3  K(  C1C4  + 
+  7  (c32  +  2C2C4  )  ;?C3C4  +  ?  C42]} 


db0  b0  da0  /  H  b0\ 

0  0 


(VII-41) 


(VII-42) 


By  using  equations  (VII-13)  through  (VII-33)  equations  (VII -40),  (VII-41) 
and  (VII-42)  may  be  reduced  to  a  system  of  three  ordinary  differential  equations 
with  dependent  variables  b^,  6^,  a°d  Since  sYstem  of  ordinary  dif¬ 

ferential  equations  involves  all  the  profile  coefficients  and  thicknesses  and  their 
derivatives  with  respect  to  x  in  a  very  complicated  manner,  these  equations 

will  be  rearranged  so  that  the  numerical  solution  can  be  conveniently  obtained 

« 

by  using  a  digital  computer.  In  general  equations  (VII -40),  (VII-41)  and  (VII-42) 
can  be  written  as 


/V  ^  ( 

where  the  s  are  functions  of  b^,  fi  6.  a^,  6,  and  q,  and  the  s  are 

functions  of  t>2,  6^,  0^  aQ,  0,  qs,  ,  -j^r  .  The  (^’s  are  giyen  in  the 


Appendix.  Defining 
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•Oie  initial  values  for  b^,  6 and  may  be  approximated  by  appropriate 
curve  fitting  of  the  similarity  velocity  profile,  for  die  mixing  between  a  uniform 


stream  and  a  fluid  at  rest,  evaluated  at  a  fixed  but  small  x  =  Xq.  Therefore, 
for  y  i  0,  equation  (VI  -43)  can  be  written  as 
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Since  b^x^)  and  c  (x^)  are  given  by  equations  (VIi-13)  through  (VII -31),  which 
are  functions  of  b2(x0),  6^(x0),  ^  _^(Xq).  a(/°^’  ttI,d  ^  since 

and  6(0)  are  known,  b.(x„;,  6, .  and 

2  0  (+)/°  ftp;  _  \ 

6(  _) (STq)  can  be  determined  approximately  from  the  plot  of  f^  f  J  y  J  . 

Hence  the  system  of  ordinary  differential  equations  (V1I-46),  (VII -47)  and  (VII -48) 
can  be  solved  numerically  by  using  the  Runge-Kutta  method. 

Since  the  shear  stress  distribution  along  die  base  of  the  wedge  and  the 
dividing  streamline  can  be  determined,  and  since  the  dividing  streamline  is 
assumed  straight,  the  conservation  of  angular  momentum  of  the  recirculation 
region  requires  that 


(f°-k 
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where  and  are  the  perpendicular  distance  from  the  base  cf  the  wecjge  and  the  dividing 
streamline  to  the  center  of  rotation  respectively,  which,  for  a  given  wake  angle, 
can  be  determined  from  the  solution  of  the  inviscid  core  flow.  From  equations 
(vn-7)  and  (VII-0),  the  above  equation  can  be  rewritten  as 
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of  state,  the  above  equation  reduces  to 
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and  the  equation 
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Since  =  a)r  L^/2  and  0)^  =  o)q,  the  equilibrium  rate  of  rotation  for  die  in- 
viscid  core  fluid  can  be  obtained  for  a  given  wake  angle: 
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Accordingly,  for  a  given  0q,  there  exists  an  equilibrium  u)q  which  in  turn  pro¬ 
vides  a  velocity  profile  for  the  shear  layer  at  x  =  1.  Hew .e  the  pressure  at  the 
rear  stagnation  point  can  be  obtained  from  the  local  Mach  number  at  x  =  1  and 

y  =  0.  Since  u  =  a0(l)bQ(l)ue  and  T  =  %  -  £  a^l)  bg(l)  J ,  then 
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and  the  recompression  ratio  for  the  dividing  streamline  is 
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2 .  Evaluation  of  Results  and  Comparison  with  Experimental  Data 

The  theoretical  model  studied  previously  concerns  the  laminar  near  wake 

behind  a  slender  wedge  as  Re  -»  »  .  However,  for  finite  but  relatively  large  R 

with  a  laminar  near  wake,  this  model  should  provide  reasonably  accurate  results . 

Computation  have  been  carried  out  for  the  flow  over  awedge  with  0c  =  6°  and 

=  0.399  under  the  following  upstream  conditions: 

M  =  2.61  and  3.51  ;  Re  =  2  X  105 
»  c 

For  these  given  conditions,  the  experimental  results  show  that  the  near  wake  is 
laminar  but  very  close  to  the  transition  range;  the  wake  angle  is  quite  small;  the 
lip -shock  strength  is  relatively  small;  the  recompression  region  is  not  too  long 
and  the  boundary -layer  thickness  is  still  small.  Above  all,  the  prevailing  rea¬ 
son  for  choosing  these  upstream  conditions  is  that  these  are  the  available  experi¬ 
mental  conditions  closest  to  the  theoretical  model,  and  a  comparison  can  be  made 
between  die  theoretical  and  the  experimental  results .  Based  on  the  above  given 
conditions,  it  can  be  shown  that,  for  a  flow  with  total  temperature  equal  to  1000°R, 
the  flow  properties  after  the  wedge  shock  are 

M  =  2.34,  T  =  477. 3°R  ,  p  =  0.678  x  10'4  slug/ft3 

c  c  c 

u  =  2500  ft/sec,  u  =  0.338  X  10'6  lb -sec/ft2  for  M  =  2.61 

c  c  m 

Mc  =  3.14,  T  =  336. 5°R  ,  P.  =  0.422  X  10'4  slu^ft3 

u  =  2810  ft/sec,  p  =  0.236  X  10’6  lb-sec/ft2  for  Mw  =  3.51 

c  c 

The  flow  properties  corresponding  to  the  conditions  at  the  outer  edge  of  the  shear 
layer  and  the  recirculation  region  may  be  obtained  provided  that  the  wake  angle 
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is  known.  The  flow  properties  after  the  corner  expansion  are 

M  =  3.01,  T  =  356°R,  p  =  0.323  x  10'4  slug/ft3 

u  =  2780  ft/sec,  p  =  0.25  X  10'6  lb-sec/ft2 
e  e 

Re  =  0.998  X  105 
e 


for  M 


6. 


2.61 

8.7° 


Mg  =  4.28,  T£  =  211°R,  pg  =  0.133  X  10*4  slug/ft3 

u  =  3074  ft /sec,  p  =  0.15  X  10  ^  lb-sec/ft2 

e  e 

Re  =  0.602  x  105 
e 

The  flow  properties  in  the  recirculation  region  are 

=  1000°R,  pr  =  0.71  X  10‘6  lb-sec/ft2 

and 


p  =  0.115  x 

I 

o 

slug/ft3 

for  M 

=  2.61 

r 

p  =  0.281  x 
r 

10'5 

slug/ft3 

m 

for  M 

OO 

=  3.51 

Since  the  complete  solution  for  the  inviscid  core  flow  is  given  by  equation 
(IV- 29),  the  streamlines  in  the  inviscid  core  flow  for  0q  =  10°are  drawn  .  Figure 
(4a)  shows  the  streamline  pattern  when  coordinates  are  nondimensionalized  by  the 
wake  length.  Figure  (4b)  shows  the  streamline  pattern  with  the  abscissa  normaliz- 
•d  by  the  length  of  the  wake  and  the  ordinate  normalized  by  the  base  height.  The 
flow  velocity  is  relatively  high  in  a  region  within  a  distance  equal  to  the  base 
height  from  the  base  of  the  wedge,  and  decreases  gradually  toward  the  rear  stag¬ 
nation  point .  It  can  be  shown  that  the  distance  between  the  center  of  rotation  and 
the  base  of  the  wedge  increases  as  9q  increases . 

Figure  5  through  figure  8  are  drawn  for  M  =  2.61,  Re  =  2  x  103, 


for  M  =  3.51 


0  =  10.7 

o 
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and  for  0  =  8.7°.  This  0  satifies  the  reattachment  condition  and  the  require  - 
o  o 

ment  of  conservation  of  angular  momemtum  in  the  recirculation  region. 

For  a  given  0q,  equations  (VII-1)  through  (VII-6)  can  provide  a  numerical  solution 
for  the  base  boundary  layer.  The  shear  stress  along  the  base  of  the  wedge  is 
shown  in  figure  5.  This  increases  very  rapidly  with  0  for  0  <  0.2,  because 
the  inviscid  core  flow  has  a  large  acceleration  for  r  =  1  and  0  <  0.2.  Since 
the  inviscid  core  flow  decelerates  for  r  =  land  0  >0.5,  the  shear  stress  di¬ 
minishes  rather  rapidly  and  approaches  7ero  near  0  =  0.8.  Figure  6  shows 
the  velocity  profiles  at  0  =  0,  0.15,  0.30,  0.45,  and  0.60.  The  boundary  - 
layer  thickness  grows  gradually  for  0  £  0.5  but  increases  very  fast  for  0  >  0.5. 
The  integral  method  used  appears  to  predict  separation,  which  may  occur  near 
0  =  0.8.  The  stagnation -point  velocity  profile  obtained  from  the  similarity 
solution  is  also  shown  in  order  that  a  comparison  can  be  made  between  the 
similarity  solution  and  the  solution  from  the  method  of  integral  relations  for 
0  =  0.  It  is  evident  that  the  method  of  integral  relations  gives  a  fairly  accur¬ 
ate  velocity  profile  for  U  <  0.9.  The  accuracy  will  certainly  be  improved,  if 
the  order  of  approximation  K  in  (V)  is  increased. 

When  0q  is  given,  the  ' itial  velocity  profile  for  the  shear  layer  can  be 
determined.  Therefore  the  numerical  solution  can  be  obtained  for  the  shear 
layer  equations  (VII -46)  through  (VII -48).  Figure  7  shows  the  shear  stress  dis¬ 
tribution  along  the  dividing  streamline .  Initially  the  shear  stress  is  very  large, 
because  the  inviscid  core  flow  velocity  is  very  small,  at  0  =  1  and  for  r  -•  1, 
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and  the  shear  sublayer  is  extremely  thin.  Since  the  inviscid  core  flow  has  large 

acceleration  for  0  =  1  and  r  >  0.85  (or  —  <  0.15),  and  since  the  thickness 

o 

of  the  shear  sublayer  grows  fairly  rapidly  for  r  >  0.85,  the  shear  stress  along 
the  dividing  streamline  declines  rather  quickly.  The  shear  stress  decreases 
slowly  for  r  <  0.85,  because  the  shear  sublayer  grows  rather  slowly  and  the 
velocity  of  the  inviscid  core  flow  decreases  gradually.  Figure  8  gives  the  velo¬ 
city  profiles  for  the  shear -layer  at  x  =  0.05,  0.2,  0.45,  0.7,  and  1.0.  For 
x  small  the  change  in  velocity  profile  is  large  near  the  dividing  streamline.  As 
x  increases  this  change  in  velocity  profile  spreads  out  to  a  larger  region  on  each 
side  of  the  dividing  streamline.  The  velocity  profile  for  the  portion  of  the  shear 
layer  with  high  kinetic  energy  changes  fairly  slowly  throughout  the  whol?  near- 
vake  region. 

Figure  9  shows  the  curves  used  to  obtain  the  wake  angle,  which  in  turn 
determines  the  base  pressure.  The  solid  curve  is  obtained  from  equation.(VII  54), 
which  furnishes  the  recompression  ratio  along  the  dividing  streamline  for  M  = 
2.61  and  for  an  equilibrium  rate  of  rotation  for  the  inviscid  core  flow.  When  a 
uniform  supersonic  flow  is  turned  through  a  given  angle  from  a  specified  initial 
static  pressure,  the  final  value  of  the  pressure  depends  on  the  initial  value  of 
the  Mach  number,  and  the  final  pressure  has  a  minimum  value  for  a  particular 
initial  Mach  number.  Hence  a  minimum  compression  ratio  can  be  determined 
and, is  shown  by  the  dotted  line.  The  intersection  between  die  solid  and  the 
49tted  curves  provides  an  appropriate  wake  angle  and  recompression  ratio, 
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which  in  turn  gives  the  base  pressure.  This  intersection  point  satifies  the  re¬ 
attachment  condition.  Since  figure  9  gives  0q  =  8.7°,  and  P/P^  =  1-52,  the 
corresponuing  equilibrium  rate  of  rotation  for  inviscid  core  flow  is  ou^  = 
-60000  rad/sec.,  and  p^/p^  =  0.538.  Hence  =  1250  ft/sec.,  Re^  = 

3 

0.844  X  10  ,  =  0.806.  The  experimental  results  given  by  Hama  for 

Re  =  2  X  10**  and  M  =  2.61  show  that  0  &  7.5°;  p./p  -  0.56and 
the  lip-shock  strength  is  approximately  equal  to  1.25.  Consequently,  the 
computed  and  the  experimental  results  agree  fairly  well  for  this  case. 

For  Re^  =  2  X  10**  and  =  3.51,  the  theoretical  model  predicts 

that  0  3  10.7°,  p./p  =  0.36,  =  1190  ft/sec.,  and  Re  =  0.196X  103 

and  =  0.78-  The  experimental  results  given  by  Hama  show  that 

0Q  sr  8.5°,  Pjj/1^  =  0.45  and  the  lip -shock  strength  is  approximately  equal 
to  1.4.  The  predicted  results  are  less  accurate  for  this  case. 


VIII 


CONCLUSIONS 

A  theoretical  description  has  been  presented  for  the  flow  in  the  laminar 
near  wake  of  a  slender  wedge  at  supersonic  speed  and  high  Reynolds  number. 

This  study  consists  of  separate  analyses  for  the  various  near -wake  regions  hav¬ 
ing  different  physical  character.  The  solutions  for  these  regions  must  be  match¬ 
ed  through  the  condition  of  the  conservation  of  angular  momentum  for  the  recir¬ 
culation  region  and  the  proposed  reattachment  condition,  so  that  the  near -wake 
flow  is  uniquely  determined. 

As  the  Reynolds  number  approaches  infinity  as  a  limit,  the  thickness  of 
the  laminar  viscous  layers  enclosing  the  recirculating  flow  is  of  order  Ree 
and  the  flow  in  the  interior  of  the  recirculation  region  may  be  considered  as  in- 
viscid  and  nonconducting.  The  typical  velocity  of  the  recirculating  flow  must  be 
of  the  same  order  as  that  of  the  external  flow,  although  numerically  it  is  con¬ 
siderably  smaller.  The  Reynolds  number  Ref  for  the  recirculation  region  dif¬ 
fers  from  Reg  through  factors  involving  powers  of  Mg.  .herefore,  in  the  limit 
Re^-.  «  with  fixed,  Ref  will  also  approach  infinity.  Since  is  fairly 
small,  in  the  zeroth -order  approximation  the  recirculating  flow  may  be  re¬ 
garded  as  incompressible.  The  inviscid  incompressible  recirculating  flow 
has  constant  vorticity  and  constant  temperature.  The  effect  of  compressibility 
is  estimated  to  be  rather  small . 

The  inviscid  incompressible  recirculating  flow  in  the  right  triangular 
region  is  characterized  by  the  small  wake  angle  0q.  Hence  the  stream  function 
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tp  may  be  expanded  in  a  power  series  in  0q,  and  an  approximate  solution  may  be 

obtained  analytically.  The  zeroth-order  solution  for  i p  is  symmetric  about  = 

the  first-order  solution  for  ip  does  not  have  this  symmetry.  It  is  not  possible 

to  obtain  a  second-or  higher -order  solution,  because  the  representation  for 

one  of  the  boundary  conditions  is  not  a  uniformly  convergent  series .  The  flow 

velocity  is  higfr  in  a  region  within  a  distance  of  the  order  of  the  base  height 

from  the  base  of  the  wedge .  The  maximum  flow  velocity  is  approximately 

0.85  q^  =  0.85  («r  L^/2),  and  the  distance  between  the  base  of  the  wedge 

and  the  center  of  rotation  increases  as  0  increases  .  It  should  be  mentioned 

o 

that  the  formulation  in  terms  of  a  right  triangular  region  assumes  that  sec¬ 
ondary  eddies  either  do  not  occur  or  have  negligible  effect. 

In  general  the  shear  stress  along  the  base  is  large  for  0  <  0.5,  be¬ 
cause  the  flow  is  accelerating.  The  shear  stress  declines  rather  rapidly  fcr 
0  >  0.5  due  to  the  adverse  pressure  gradient.  The  boundary-layer  thickness 
grows  gradually  for  0  <  0.5  and  more  rapidly  for  0  >  0.5.  Although  the  in¬ 
tegral  method  predicts  boundary -layer  separation  near  0  =  0.8,  it  is  believed 
that,  due  to  the  suction  effect  of  the  shear  layer  near  the  trailing  edge  of  the 
wedge,  the  separation  will  be  delayed. 

The  development  of  the  shear  layer  is  investigated  by  including  the 
effects  of  the  finite  initial  thickness,  the  trailing-edge  expansion,  the  nonzero 
innet-edge  velocity,  and  the  equalization  of  the  pressure  across  the  shear  layer 
up  to  0(ReeV  However,  it  is  not  necessary  to  consider  the  details  of  the  cor- 
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ner  expansion,  because  only  a  small  distance  is  required  for  the  completion  of 
the  expansion,  and  because  the  lip-shock  strength  is  rather  small  (pressure 
ratio  1 .25  for  Mw  =  2.61  and  1 .40  for  =  3.51 ).  The  portion  of  the  vel¬ 
ocity  profile  with  high  kinetic  energy  remains  nearly  unchanged  throughout  the 
near -wake  region,  but  die  portion  of  the  velocity  profile  near  the  dividing 
streamline  is  changed  greatly.  These  results  suggest  that  it  is  indeed  appro¬ 
priate  to  divide  the  shear  layer  into  an  outer  shear  layer  and  a  shear  sub¬ 
layer,  and  that  the  changes  in  the  outer  shear  layer  actually  could  have  been 
neglected.  Initially  the  thickness  of  the  shear  sublayer  grows  rapidly,  and 
the  shear  stress  along  the  dividing  streamline  declines  rapidly. 

The  proposed  reattachment  condition  is  based  on  the  minimum  pres¬ 
sure  increase  attainable  at  the  rear  stagnation  point  from  the  turning  of  a 
streamline  in  the  shear  layer  to  a  direction  parallel  to  the  wake  center  line. 
This  condition  and  the  requirement  that  angular  momentum  be  conserved  in 
the  recirculation  region  are  used  to  obtain  a  unique  solution  for  the  base  pres¬ 
sure  and  hence  for  the  near -wake  flow.  For  M  =  2.61  and  Re  =  2  X  10  , 

oo  c 

the  predicted  values  for  0q  and  P^/P,*,  are  fairly  close  to  the  experimental 

2  5 

results".  However,  for  M  =  3.51  and  Re  =  2  X  10  ,  the  predicted  val- 

oo  c 

ues  are  less  satisfactory,  because  of  the  smaller  value  of  Re^. 

The  deviation  of  the  predicted  values  for  M  =2.61  from  the  experi- 
mental  results  is  probably  caused  mainly  by  the  inadequacy  of  the  proposed  re- 
attachment  condition.  The  experimental  results  show  that  the  flow  is  only 


-Im¬ 


partially  recompressed  at  reattachment.  The  theoretical  reattachment  condition 

9 

proposed  by  Chapman  assumes  that  the  recompression  is  completed,  and  is 

bound  to  result  in  disagreement  between  the  predicted  and  the  measured  values 

of  the  base  pressure.  The  present  reattachment  condition  is  intended  to  remedy 

tnis  inadequacy.  However,  since  the  predicted  base  pressure  is  too  low  it 

seems  that  this  reattachment  condition  does  not  provide  sufficient  recompression. 

By  requiring  tha  ‘  a  different  streamline  be  turned  to  a  direction  parallel  to  the 

wake  center  line,  one  can  obtain  values  of  fl  and  p,  /  p  which  are  significant- 

o  b  * 

ly  closer  to  the  experimental  results . 

It  is  known  that  the  separation  along  the  base  occurs  somewhere  below 
the  trailing  edge,  that  the  dividing  streamline  is  slightly  curved,  and  that  the 
center  of  rotation  is  displaced  because  of  the  base  boundary  layer.  The  errors 
which  result  from  neglecting  these  effects  will  influence  the  geometry  of  the  in- 
viscid  recirculating  flow.  The  corresponding  changes  in  0q  and  pfa  arise  pri¬ 
marily  from  a  change  in  the  length  of  the  near  wake,  and  can  be  shown  to  be 
quite  small,  of  the  order  of  a  few  per  cent. 

For  =  2.61  and  Re^  =  2  x  1CT\  the  thickness  of  the  base  boundary 
layer  at  a  distance  of  half  the  base  height  from  the  wake  center  line  is  approxi¬ 
mately  0.07  "  ,  and  the  thickness  of  the  lower  shear  sublayer  at  x  =  £  is  approx¬ 
imately  0.05  ".  Since  the  base  height  is  0.5  ",  these  values  indicate  that  there 
exists  a  quite  large  region  of  inviscid  recirculating  flow .  Even  for  M  =3.51 

w  OO 

and  Rec  =  2  x  10^  a  significant  portion  of  the  recirculating  flow  may  be  des- 
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cribed  as  approximately  inviscid.  Therefore,  for  a  relatively  large  Reynolds 
number  the  present  model  does  give  a  fairly  good  representation  of  the  actual 
flow.  For  large  Ms  ,  however,  the  trailing-edge  expansion  region  and  the  re¬ 
compression  region  are  no  longer  small,  the  lip  shock  is  stronger,  the  wake 

Reynolds  number  is  smaller,  and  the  wake  angle  is  larger.  These  features 

19 

are  all  incorporated  in  the  analysis  carried  out  by  Weiss. 

In  order  to  understand  the  near -wake  flow  more  thoroughly  it  would  be 
necessary  to  study  the  recompression  and  the  trailing-edge  expansion  regions. 
As  mentioned  before,  the  present  reatcachment  condition  does  not  give  suffi¬ 
cient  pressure  recovery.  To  obtain  a  more  accurate  reattachment  condition 
and  to  verify  the  assumption  of  isantropic  recompression  within  a  narrow 
region,  if  would  be  necessary  to  study  the  recompression  region  in  more  de¬ 
tail.  The  separation  processes,  the  behavior  of  the  base  boundary  layer 
approaching  the  separation  point,  and  the  formation  of  the  lip  shock  could 
be  better  understood  through  the  investigation  of  the  trailing-edge  expansion 
region.  The  dependence  of  p^  and  0q  on  Reynolds  number  Joes  not  appear 
in  the  current  approximation.  A  first  approximation  for  this  dependence 
could  presumably  be  derived  after  studying  the  recompression  region  and 
the  trailing-edge  expansion  region . 


Trailing -edge  Expansion  Fan 
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Figure  2:  Proposed  Model  for  the  Near  Wake  of  Slender  Wedge  in  Supersonic  Flow 


Streamline  Pattern  with  Stretched  Vertical  Coordinate 
Streamlines  in  the  Inviscid  Core  Flow 


Compression  Ratio  P/  P 


Figure  9:  Determination  of  a  Unique  Compression  Ratio  for  the 
Dividing  Streamline 
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